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“ Surface  current  and  charge  densities  induced  on  aircraft  fuselage  by  an 
electromagnetic  plane  wave  of  arbitrary  incidence  at  high  frequencies  is  the 
object  of  this  research.  The  primary  goal  of  this  investigation  is  to  develop 
a theoretical  solution  for  predicting  the  induced  surface  and  charge  densities 
on  aircraft  fuselage  In  an  accurate  and  efficient  manner.  Since  it  is  a study 
of  general -type  aircraft,  the  aircraft  is  modeled  in  its  most  basic  form.  The 
fuselage  is  assumed  to  be  an  infinitely  long,  perfectly  conducting  -)  (over) 
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elliptic  cylinder  in  its  cross-section  and  a composite  elliptic  cylinder 
in  its  elevation  profile.  The  wing,  cockpit,  stabilizers  (horizontal  and 
vertical)  and  landing  gear  are  modeled  by  ^n"  sided  bent  or  flat  plates 
which  can  be  arbitrarily  attached  to  the  fuselage. 

The  solution  developed  in  this  study  utilizes  two  elliptic  cylinders, 
namely,  the  roll  plane  and  elevation  plane  models  to  approximate  the  princi- 
pal surface  profile  (longitudinal  and  transverse)  at  the  observation  location. 
With  the  belt  concept  and  the  aid  of  appropriate  coordinate  system  trans- 
formations the  solution  can  be  used  to  predict  the  surface  current  and  charge 
densities  induced  on  the  fuselage  in  an  accurate  and  efficient  manner.  The 
radiation  patterns  due  to  infinitesimal  monopole  and  slot  antennas  on  aircraft 
are  related  to  the  surface  current  and  charge  densities  induced  on  the  antenna 
location  via  the  reciprocity  theorem. 
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SECTION  I 


INTRODUCTION 

The  object  of  the  research  is  to  calculate  the  surface  current 
and  charge  density  induced  on  aircraft  by  an  incident  electromagnetic 
plane  wave.  The  primary  goal  of  this  investigation  is  to  develop  a 
theoretical  solution  for  predicting  the  induced  surface  current  and 
charge  densities  on  the  aircraft  fuselage  in  an  accurate  and  efficient 
manner.  Since  it  is  a study  of  general -type  aircraft,  the  aircraft  is 
modeled  in  its  most  basic  form.  The  fuselage  is  assumed  to  be  an  in- 
finitely long  perfectly-conducting  elliptic  cylinder  in  its  cross-section 
and  a composite  elliptic  cylinder  in  its  elevation  profile.  The  wing, 
cockpit  and  stabilizers  (horizontal  and  vertical)  are  modeled  by  "n"- 
sided  bent  or  flat  plates  which  can  be  arbitrarily  attached  to  the 
fuselage.  The  cylinder  solution  employed  in  the  analysis  is  quite 
accurate  and  useful  in  that  reasonable  approximations  are  made  in  ob- 
taining the  necessary  diffraction  terms. 

The  surface  current  and  charge  density  induced  on  the  aircraft 
fuselage  with  plane-wave  incidence  are  related  to  the  radiation  patterns 
due  to  fuselage-mounted,  infinitesimal  monopole  and  slot  antennas.  The 

relation  between  them  will  be  discussed  in  the  following.  Let  (E°,E^), 

0,0,  m m' 

(E  .E;)  and  (E°,E;)  be  the  radiation  patterns  due  to  an  infinitesimal 

monopole,  axial  slot  and  circumferential  slot  antenna  respectively. 

Let_  (p^1 -.p^1 ) be  the  induced  surface  charge  density,  and  (J^1  ) , 

(J^  .J^1)  be  the  induced  surface  current  density  due  to  e-  and  <j>- 

incident  plane  waves  with  unit  electric  field  intensity.  From  the 
reciprocity  theorem,  the  following  relations  exist: 
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Note  that  the  fuselage  axis  is  assumed  to  be  the  z-axis  of  the 
conventional  spherical  coordinate  system  where  o (theta)  is  the  angle 
measured  from  the  z-axis  and  <p  (phi)  is  the  angle  measured  from  the 
x-y.  plane.  A e-incident  plane  wave  means  tbe  electric  field  intensity 
associated  with  the  incident  plane  wave  is  e-polarized.  Similarly,  a 
^-incident  plane  wave  means  the  electric  field  intensity  is  ^-polarized. 

Since  extensive  experimental  scale-model  measurements  for  flush- 
mounted  and  blade  antennas  on  aircraft  are  available  for  comparison,  the 
analysis  for  the  induced  surface  current  and  charge  density  is  carried 
out  in  terms  of  antenna  problems. 

Based  on  the  two  principal  plane  (roll  and  elevation)  model  analyses 
from  previous  work  [1]  and  extensive  work  on  prolate  spheroids  [2],  a 
numerical  solution  for  the  volumetric  pattern  of  fuselage  mounted 
antennas  is  developed.  The  solution  present  here  not  only  predicts  the 
volumetric  pattern  accurately,  it  also  reduces  computer  time  and  improves 
computational  efficiency  considerably  over  that  of  a previous  surface  of 
revolution  approach  [2].  The  validity  and  capability  of  this  new  solu- 
tion are  illustrated  through  comparison  with  extensive  scale-model 
measurements  performed  at  NASA  (Hampton,  Va.). 

Among  the  first  solutions  used  to  compute  on-aircraft  antenna  pat- 
terns were  the  modal  solutions  for  infinitely  long  circular  [4,5]  and 
elliptical  [6]  cylinders.  These  solutions  modeled  the  fuselage  by  a 
cylinder  whose  elliptical  cross-section  approximated  the  fuselage  cross 
section  at  the  antenna  location.  Arbitrary  antennas  were  considered  in 
these  studies  in  which  the  antenna  was  mounted  either  on  or  above  the 
fuselage.  However,  these  solutions  were  not  always  adequate  in  that  the 
effects  of  various  scatterers  such  as  wings  and  stabilizers  were  ignored. 

A more  recent  approach  for  solving  antenna  pattern  and  impedance 
problems  is  the  integral  equation  method  via  moment  methods.  By  enforc- 
ing the  boundary  conditions  on  aircraft  structure,  the  surface  currents 
and  the  resulting  scattered  fields  can  be  found.  One  of  the  first  moment 
solutions  applied  to  aircraft  problems  was  the  wire  grid  technique,  de- 
veloped by  Richmond  [7],  using  a point  matching  scheme  [8].  This  solution 
required  the  determination  of  approximately  100  unknown  currents  per 
square  wavelength  in  order  that  the  wire  grid  adequately  model  a per- 
fectly conducting  surface.  A more  sophisticated  approach  has  been  developed 
by  Richmond  [9]  in  which  the  reaction  technique  is  used  to  solve  for  the 
unknown  currents.  Yet  this  solution  still  requires  a wire  grid  model  of 
the  aircraft  with  approximately  100  unknown  currents  per  square  wave- 
length. 

Another  approach  is  the  surface  current  model  method  [10]  in  which  the 
surface  of  a conducting  body  is  divided  into  patches  with  each  patch 
having  two  orthogonal  unknown  complex  currents.  Using  this  approach,  the 
unknown  currents  have  been  reduced  to  about  20  per  square  wavelength. 

Thus,  one  is  able  to  consider  a much  larger  surface  using  these  surface 
patch  solutions.  However,  all  of  these  solutions  are  restricted  to  lower 
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frequencies,  based  on  the  fundamental  limitation  on  the  size  of  matrices 
which  modern  computers  can  solve  without  excessive  loss  of  accuracy. 

Another  approach  that  has  found  great  success  at  analyzing  aircraft 
antennas  is  the  Geometrical  Theory  of  Diffraction  (GTD).  GTD  is  basically 
a high  frequency  solution  which  is  divided  into  two  basic  problems;  these 
being  wedge  diffraction  and  curved  surface  diffraction.  The  wedge  dif- 
fraction solution  has  been  applied  to  determine  the  radiation  patterns 
of  such  basic  antennas  as  parallel  plate  antennas  [11,12,13],  parallel 
plate  arrays  [14,15],  horn  antennas  [16,17],  parabolic  reflectors  [13,19], 
and  rectangular  waveguide  antennas  [20].  Both  of  these  diffraction 
solutions  have  been  applied  in  computing  the  patterns  of  antennas  mounted 
on  cylinders  [21,22,23],  rockets  [24],  wings  [25,26],  and  aircraft  [1,2, 
27-32].  Using  this  approach,  one  applies  a ray  optics  technique  to  de- 
termine components  of  the  field  incident  on  the  various  scatterers.  Com- 
ponents of  the  diffracted  field  are  found  using  the  GTD  solutions  in  terms 
of  rays  which  are  summed  with  the  geometrical  optics  terms  in  the  far 
field.  The  rays  from  a given  scatterer  tend  to  interact  with  the  other 
structures  causing  various  higher  order  terms.  In  this  way  one  can  trace 
out  the  various  possible  combinations  of  rays  that  interact  between 
scatterers  and  determine  and  include  only  the  dominant  terms.  Thus,  one 
need  only  be  concerned  with  the  important  scattering  components  and 
neglect  all  other  higher-order  terms.  This  makes  the  GTD  approach  ideal 
for  a general  high  frequency  study  of  aircraft  antennas  in  that  only  the 
most  basic  structural  features  of  the  aircraft  need  to  be  modeled. 

Since  GTD  is  essentially  a high  frequency  solution,  the  lower  fre- 
quency limit  of  this  solution  is  dictated  by  the  spacings  between  the 
various  scattering  centers  in  that  they  should  be  at  least  a wavelength 
apart.  In  some  cases  even  this  requirement  can  be  relaxed.  Under  this 
restriction,  the  low  frequency  limit  is  typically  around  100  MHz.  The 
upper  frequency  limit  is  dependent  on  how  well  the  theoretical  model 
simulates  the  important  details  of  the  actual  structure. 

The  basic  approach  applied  in  the  present  study  is  based  on  previous 
work  [1]  which  demonstrated  the  capability  of  the  numerical  solutions  to 
predict  the  radiation  patterns  of  fuselage  mounted  antennas  in  an  ef- 
ficient and  economical  way.  If  the  volumetric  patterns  were  found 
directly  by  analyzing  rays  on  complex  three-dimensonal  surfaces  as  done 
previously  in  Reference  [2],  the  resulting  numerical  solution  would  be 
very  complex,  time-consuming,  and  uneconomical.  Nevertheless,  if  certain 
assumptions  can  be  made,  the  approach  undertaken  previously  can  be  used  to 
overcome  these  difficulties  and  simplify  the  problem  a great  deal. 

First,  it  has  been  shown  by  comparison  with  numerous  scale  model 
measurements  that  the  roll  plane  model  (an  infinite  elliptic  cylinder 
with  flat  plate  wings)  can  be  extended  to  cover  almost  the  complete 
volumetric  pattern  except  for  two  conical  sectors  (fore  and  aft)  [30]. 

The  limitations  of  the  roll  plane  model  are  due  to  the  finite  length  of 
fuselage.  Yet,  the  problem  of  finite  length  fuselage  has  been  solved, 
previously,  in  an  elevation  plane  analysis  (an  infinite  composite 
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elliptic  cylinder).  Furthermore,  based  on  previous  three-dimensional 
studies  of  geodesic  rays  which  contribute  to  the  radiation  pattern  of  an 
antenna  on  various  prolate  spheroids,  one  is  able  to  combine  the  analysis 
of  these  two  models  to  give  the  complete  pattern. 

In  addition,  the  cockpi t/radome  section  and  vertical  stabilizer 
previously  ignored  are  taken  into  account  in  the  present  study.  The 
cockpi t/radome  section  and  stabilizers  are  approximated  by  flat  or  bent 
plates  which  can  be  arbitrarily  attached  to  the  fuselage.  This  is  an 
improvement  over  the  previous  roll  plane  model  [1].  Using  this  approach, 
the  complete  volumetric  pattern  can  be  obtained  with  a model  oOnSistiny 
of  a composite  elliptic  cylinder  to  which  are  attached  flat  or  bent 
plates.  As  a result  of  this  simplified  model,  the  solution  is  very 
efficient  and  requires  little  computer  storage. 

The  basic  theoretical  background  on  the  geometrical  theory  of 
diffraction  (GTD)  is  presented  in  Section  II.  Both  wedge  diffraction 
and  curved  surface  diffraction  are  discussed.  The  infinitely  long 
elliptic  cylinders  (curved  surfaces)  is  also  presented  in  that  section. 
This  formulation  is  essential  for  fuselage  mounted  antennas.  Section  II 
also  includes  the  near  field  scattering  by  a finite  bent  plate  which  is 
necessary  for  the  determination  of  the  scattering  effects  of  secondary 
contributors  such  as  wings,  cockpit,  and  stabilizers. 

Section  III  describes  the  mathematical  model  of  the  aircraft. 

This  section  starts  with  a review  of  the  two  principal  plane  (roll  anc 
elevation)  analyses  and  their  application  to  the  radiation  patterns. 

It,  then,  explains  the  way  these  two  analyses  are  combined  into  a com- 
plete solution  for  the  three  dimensional  aircraft  model.  The  computer 
simulation  method  to  approximate  the  cross  section  and  elevation  profile 
of  an  aircraft  fuselage  is  also  discussed. 

Section  IV  presents  the  numerical  results  for  the  antenna  patterns 
using  these  newly  developed  solutions.  Volumetric  patterns  of  both  cal- 
culated and  the  experimental  results  are  presented  to  demonstrate  the 
validity  and  capability  of  this  solution.  Section  V presents  the  calcu- 
lated surface  current  and  charge  densities  induced  on  various  aircraft 
models.  A complete  volumetric  pattern  for  the  charge  density  induced  on 
a Boing  737  is  also  presented  in  this  section.  Finally,  a summary  of 
the  present  study  and  a discussion  of  future  topics  are  presented  in 
Section  VI . 
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SECTION  II 


THEORETICAL  BACKGROUND 


As  stated  in  the  introduction,  the  Geometrical  Theory  of 
Diffraction  (GTD)  is  ideal  for  a general  high  frequency  study  of 
fuselage  mounted  aircraft  antennas  and  induced  surface  densities  of 
current  and  charge.  This  is  particularly  true  when  the  scattering 
object  is  large  in  terms  of  the  wavelength  such  as  for  electrically 
large  aircraft.  An  aircraft  shape  is,  in  general,  quite  cumplex  con- 
sisting of  many  complicated  scattering  structures.  To  be  able  to  obtain 
an  accurate  radiation  pattern,  one  must  take  these  various  scattering 
structures  into  account. 

Based  on  past  performance,  the  GTD  has  proven  itself  well  suited 
to  this  type  of  analytical  study.  Not  only  does  this  approach  fit 
nicely  into  a ray  optics  format,  but  it  also  provides  a means  of 
analyzing  the  effects  of  three-dimensional  structures  and  identifying 
the  significant  contributions  in  the  resulting  antenna  pattern.  Con- 
sequently, the  GTD  is  employed  in  this  study  to  analyze  the  volumetric 
patterns  of  fuselage-mounted  antennas. 

The  Geometrical  Theory  of  Diffraction  was  introduced  by  Keller  [33] 
as  an  extension  of  geometrical  optics  to  include  diffracted  fields  in  the 
high  frequency  solution.  The  theory  is  based  on  the  following  postulates: 

(1)  The  diffracted  field  propagates  along  rays  which  are 
determined  by  a generalization  of  Fermat’s  Principle  to 
include  points  on  edges,  vertices,  and  smooth  surfaces  in 
the  ray  trajectory. 

(2)  Diffraction  like  reflection  and  transmission  is  a local 
phenomenon  at  high  frequencies,  i.e.,  it  depends  only  on 
the  nature  of  the  boundary  surface  and  the  incident  field 
in  the  immediate  neighborhood  of  the  point  of  diffraction. 

(3)  The  diffracted  wave  propagates  along  its  ray  so  that 

(a)  power  is  conserved  in  a tube  (or  strip  of  rays), 

(b)  the  phase  delay  along  the  ray  path  equals  the 
product  of  the  wave  number  of  the  medium  and 
the  distance  s. 
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Using  these  postulates,  one  can  express  the  diffracted  field 
in  the  same  form  as  a geometrical  optics  field  with  some  coef- 
ficient of  proportionality  to  the  incident  field  at  the  point  of 
diffraction.  The  coefficient  is  determined  from  a canonical 
problem  and  is  referred  to  as  a diffraction  coefficient. 

For  practical  purposes,  the  GTD  can  be  divided  into  two 
categories:  (1)  wedge  diffraction  theory  - to  treat  diffraction 
by  edges  and  (2)  creeping  wave  theory  - to  treat  diffraction  by 
curved  surfaces.  The  GTD  has  evolved  considerably  from  its 

pi  ral  ■frv  v*m  in  Dofo^arirnr  TQ/ll  d n H Pnnr  A/iiinnf  1 «/ 

w*  ' -5  ■ ' ° H..W.  N.M1.V.J  L ^ ' J V*"'-*  L ^ J • vUius.vjuLiM,iy» 

descriptions  of  the  wedge  diffraction  problem  and  the  diffraction 
by  a curved  surface  are  presented  in  the  following  paragraphs. 


1.  DIFFRACTION  BY  A WEDGE 

An  asymptotic  solution  for  the  diffraction  from  a con- 
cucting  wedge  was  first  solved  by  Sommerfeld  [36].  Originally, 
plane  wave  diffraction  coefficients  as  presented  by  Keller  were 
used  as  the  sole  GTD  solution;  however,  as  shown  in  Reference 
[37]  diffraction  of  cylindrical  waves  is  necessary  in  the  treat- 
ment of  antennas. 

Pauli  [33]  introduced  the  Vg  function  as  a practical 
formulation  to  the  solution  for  a finite  angle  conducting 
v.ecge.  Recently,  however,  Hutchins  and  Kouyoumjian  [39,40] 
presented  a formulation  for  the  diffracted  field  (Vg) , which 
significantly  improves  the  accuracy  over  that  obtained  from 
^ajli's  form.  This  improved  diffraction  solution  [39,40]  Vg ( L ,<i> , n ) 
provides  superior  results  in  the  transition  regions  (near  the 
incident  and  reflected  shadow  boundaries).  It  can  be  written 
in  the  form 


VB(L,$,n)  = I_7t(L,<j»  ,n)  + 1^(1.,$, n)  (1) 

where 

-j  (kL  + tt/4)  / \ 

I±n(L,$,n)  = — —=  /a  cot(^J  x (2) 

* k I f 00  2 

d j e dx  + [higher-order  terms] 

jiU 
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and  where  the  higher-order  terms  are  negligible  for  large  kL 
and  with  n defined  by  the  wedge  angle  [WA  = (2-n)ir],  Also, 
a = 1 + cos  U-2nnN)  and  N is  a positive  or  negative  integer 
or  zero,  whichever  most  nearly  satisfies  the  equations 


2nnN-i{i  = -it 


2n7iN-4i  = +ii 


for  I 

-IT 

(3a) 

for  I 

+ TT 

(3b) 

Tne  variables  L and  $ are  defined  later. 

The  three-dimensional  wedge  diffraction  problem  is  depicted 
ir.  Fig.  1.  A source  whose  radiated  E field  is  given  by  Tn(s)  is 
located  at  point  s'.  It  can  be  an  arbitrary  electric  or  magnetic 
source  causing  a plane,  cylindrical,  conical,  or  spherical  wave 
irciaer.ce  on  the  wedge.  The  diffracted  vector  field  at  observa- 
tion point  s can  be  written  in  terms  of  a dyadic  diffraction 
coefficient.  Kouyoumjian  and  Pathak  [41,42]  have  given  a more 
rigorous  basis  for  the  GTO  formulation  and  have  shown  that  the 
diffracted  fields  may  be  written  compactly  if  they  are  in  terms 
of  a ray-fixed  coordinate  system.  The  ray-fixed  coordinate 
system  is  centered  at  the  point  of  diffraction  Qe  (or  points  of 
c,  ffraction  in  the  case  of  plane  wave  incidence).  Qe  is  a unique 
point  (or  points)  for  a given  source  and  observation  point.  The 
incident  ray  diffracts  as  a cone  of  rays  such  that  the  half  cone 
angle  e0  = S'0>  the  angle  which  the  incident  ray  makes  with  the 
edge  . 

The  relationship  between  the  orthogonal  unit  vectors  associ- 
ated with  the  ray-fixed  coordinate  system  (s', 6'  s,fc  ,$)  are 

given  by  o o 

*1  = :s' 

1 = S'q  xy  (4) 

s - Sq  * ^ 


wrere  I is  the  incident  direction  unit  vector,  and  s is  the 
diffraction  direction  unit  vector.  The  diffracted  field  is, 
now,  given  by 

^(s)  = Ei(QE)  • irE(s,i)  A(s)  e’jks  (5) 

The  quantity  A(s)  is  defined  in  Equation  (13.)  and  Qe  is  the  diffraction 
point  located  on  the  edge. 
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T 


OBSERVATION  POINT 


Figure  1.  Geometry  for  three-dimensional  wedge  diffraction 
problem. 
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Where 


De(s,I)  = -Bq  BqOs  - Dh 


Equation  (6)  is  a dyadic  diffraction  coefficient  as  given  in 
Reference  [41,42].  The  quantities  Ds  and  are  the  scalar  diffrac- 
tion coefficients  for  the  soft  (Dirichlet)  boundary  condition  at 
the  surface  of  tne  wedge  and  the  hard  (Neumann)  boundary  condition, 
respectively.  These  diffraction  coefficients  Ds  and  are 
related  to  the  Vg  function  described  earlier  by  a constant  as 


J TejkL 


Vr.  + 

sin  Bq  B 


where  L will  be  defined  later  when  applied  to  straight  edge. 

For  our  purpose,  it  is  more  convenient  to  write  the  diffracted 
field  in  terms  of  the  Vg  function  in  Eq.  (7)  as* 


E,  (s) 


E?(s). 


o -vj]  Le;(qe) 


e:,(qe) 


A,s,  .-1“ 


sin  6 


where 


Vg  = Vg(L,4>"  ,n)  + VB(L,*+,n) 


The  minus  sign  (Vg)  applies  for  the  E-field  vector  parallel  to  the 
edge  with  boundary  condition 


^ wedge ^ = 0 


*If  an  edge  fixed  coordinate  system  is  used,  Eq.  (8)  takes  the  form 
of  a 3 x 3 matrix. 


I 


In  the  far  field  (p>>p'),this  becomes 


EjJ  (p  ,<f) 

'V, 

-VB  0 ' 

Efi(p ' »4» ' ) 

_E^(p,*)_ 

- 0 -V+B- 

JT 


jko'  1 


•jkp 


J*  (16) 


e- jkp 

Putting  this  in  ray  form  and  factoring  out  -w  > one  obtains 

J p 


R?,(4>) 

-VB  0 

rJu-) 

rJu) 

o -v+B 

RjU') 

(17) 


The  ray  form  used  here  is  given  by 

-jkp 


El  ( p , <t> ) = R(f>) 


JT 

Thus,  R($)  is  related  to  the  far  field  pattern  function. 

For  the  three-dimensional  wedge  problem,  where  there  is 
spherical  wave  incidence,  Eq.  (8)  reduces  to  give 


(18) 


Ej|(s, 

-VB 

0 

'U 

0 

-v+ 

VBJ 

Elfs'.s;,*') 

2 

jks's  sin  B 


s'  s + s' 

s+s e 


o e-Jks 


(19) 


In  the  far  field  (s>>s‘),  we  have 


Eu  ( s , Bq  »^  )" 

'V, 

‘-VB 

0 ' 

"e,|  ( s 1 , Bq  1 ) 

2 

jks'sin  bq  e-jks 
s'e  s 

-E^(s,B0»^), 

- 0 

-V+ 

VBJ 

.E^s'.b’,^] 

(20) 
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~,r  - 


Putt i nq  this  in  ray  form  and  factoring  out 


, it  is  soon  that 


p-jks 


--v-  0 

% 

I 

-e- 

o 

QQ 

- 0 -VB  - 

g-jks'  cos2Bq 

(21) 


It  is  interesting  to  note  that  in  the  principal  plane  (B0  = 90°) 
the  ray  form  of  the  three-dimensional  case  takes  on  the  same  form 
as  the  two-dimensional  problem. 

The  total  ray  value  at  the  observation  point  s is  given  by 
the  sum  of  the  geometrical  optics  terms  and  the  diffracted  terms 


RT(s)  = R^°-(s)  + R^s) 


(22) 


where 


R 1 ( s ) + R^s) 


R°- 


0. 


(s)  ^ 


R’(s) 


0 


incident  and  reflected  region  I 
incident  region  II  (23) 

shadowed  region  III 


and  R (s)  may  be  determined  from  the  image  of  the  source  term  using 
basic  geometrical  optics  tecnniques.  These  three  regions  are 
illustrated  in  Fig.  2 for  a two-dimensional  wedge  diffraction 
problem. 


2.  DIFFRACTION  BY  A CURVED  SURFACE 

When  an  incident  ray  strikes  a smooth,  curved  perfectly 
conducting  surface  at  grazing  incidence,  i.e.,  at  the  shadow 
boundary,  a part  of  its  energy  is  diffracted  into  the  shadow 
region.  To  describe  this  phenomenon,  Keller  [43]  introduced 
a second  class  of  diffracted  rays  which  is  now  well  known  as 
creeping  waves.  These  ray  paths  include  the  points  Q]  and  Qo 
which  form  a curve  on  the  diffracting  surface,  as  illustrated 
in  Fig.  3.  However,  the  actual  concept  of  creeping  waves  was 
introduced  by  Franz  and  Depperman  [44,45].  The  basic  concept 
as  presented  in  the  following  discussion  is  taken  from 
"Asymptotic  High-Frequency  Methods"  by  Kouyoumjian  [46]. 
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WAVEFRONT 

TOP  VIEW 


SHADOW 

BOUNDARY 


SIDE  VIEW 


Figure  3.  Diffraction  by  a smooth  curved  surface. 


The  diffraction  by  a smooth  curved  surface  is  shown  in 
Fig.  3 in  which  0 is  the  source  point  and  P is  the  observation 
point  in  the  shadow  region.  Applying  Fermat's  principle,  the 
line  OQ1Q2P  is  the  shortest  distance  between  0 and  P which  does 
not  penetrate  the  surface.  In  detail,  a ray  incident  on  the 
snadow  boundary  at  Q]  divides;  one  part  of  the  incident  energy 
continues  straight  on  as  predicted  by  geometrical  optics,  a 
second  part  follows  the  surface  s into  the  shadow  region  as  a 
surface  ray  shedding  diffracted  rays  tangentially  as  it 
propagates  where  ?,  ft,  and  & are  the  unit  vectors  in  the 


cirection  of  incidence,  normal  to  the  surface  s and  binormal 
to  tne  surface  (6  = t x n),  respectively.  The  incident  field 
^(QlLmay  be  resolved  into  its  normal  and  binormal  components 
(n  • E'(Qi)  and  b • E 1 ( Qi ) ) . It  is  assumed  that  these  two 


components  induce  surface  ray  fields  which  propagate  inde- 
pendently of  each  other  along  the  geodesic  arc  between  Qi  and 

q2. 
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From  Reference  [46]  the  binormal  surface  ray  field  at  Q] 
is  related  to  the  binormal  component  of  the  incident  field  at  Qi 

by 


j 4>  ( Q I ) a • 

A(Q1)  e = DS(Q1 ) b]  • E1 (Q] ) (24) 

wnere  ( Q 1 ) is  the  scalar  diffraction  coefficient  for  a soft 

surface.  The  amplitude  of  the  surface  ray  is  assumed  to  be 
governed  by  the  conservation  of  energy  between  a pair  of  adjacent 
surface  rays.  Hence,  the  amplitude  behavior  of  the  fields  is 
given  as 


A(Q2)  = A(Q1 ) 


- 

rQ  2 

- 

cx(t')  dt ' 

- 

Q1  J 

(25) 


wnere 

dr, ^ and  dn2  = the  separation  between  a pair  of  rays  at 
and  respectively. 

a(t)  = the  attenuation  constant  wich  is  a function 
of  t,  the  coordinate  along  the  surface  ray, 
because  it  depends  on  the  local  radius  of 
curvature  and  its  derivatives. 


Tne  attenuation  constant  a(t)  is  introduced  due  to  the 
tangential  sheddinq  of  rays  as  the  surface  rav  propaqates.  It 
is  seen  from  Fig.  3 that  Q?  is  a caustic  of  the  diffracted  field 
and  the  second  caustic  is  located  at  a distance  p from  Q2  Thus, 
the  binormal  component  of  the  diffracted  field  which  radiates 
from  Q2  towards  P can  be  found,  as  in  the  previous  edge  diffrac- 
tion case,  with  one  of  the  caustics  used  as  a reference  point  and 
is  given  by 


eV)  = 


W 


a(q2) 


J$(Q?) 


• it 

JsCp 


£ _ 

+ s7 


,-jks 


(26) 
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From  Eqs.  (24),  (25),  and  (26),  there  results 


b2  • eV) 


b,  • E1 (Q1 ) 


Os(Qi)  VV 


p 

s(p  + s) 


-j|k(t+s)  + 
e 


Q2 

a(t')  dt' 


(27) 


It  is  found  that  • E1 (Qi)  excites  an  infinity  of  surface 
ray  modes  each  with  its  own  diffraction  coefficient  and  attenua- 
tion constant.  Thus,  the  expression  in  Eq.  (27)  is  replaced  by 


b2  • EV)  = 


dn-i 


b1  ' £ (VjJ4  sTTTJJ 


^-3  Ck( t+s ) ] 


«Jf)  dt' 


(28) 


Equation  (23)  relates  the  diffracted  field  at  P to  the  incident 
field  at  for  the  soft  surface  boundary  condition. 

An  expression  similar  to  Eq.  (23)  is  also  obtained  for  the 
normal  component  of  the  incident  field;  in  this  case,  the  scalar 
diffraction  coefficients  and  attenuation  constants  for  the  hard 
surface  replace  those  of  the  soft  surface.  Therefore,  the  vector 
diffracted- field  at  P can  be  written  in  terms  of  the  electro- 
magnetic field  incident  at  Q1  as 


["2" 


E(P) 

in  which  v(l,2),  u(l,2)  are 

v(l ,2)  = U(1 ,2)  = 


, .0,2)  * b2b,  u(l,2)l  ■ 


Jdn 


dn. 


1 Q-jkt 


[DJD  DJ2)  e 


ajf)  df 


;-jks 

(29) 


(30) 
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with  the  subscripts  h,  s,  respectively,  added  to  Dm  and  am. 

Note  that  Q-j  and  Q2  have  been  replaced  by  1 and  2 for  the  sake 
of  brevity.  Finding  dn] , dn2>  and  p is  simply  a matter  of 
differential  geometry  involving  the  rays  and  surface;  this  is 
discussed  at  length  in  Levy  and  Keller  [43],  The  generalized 
diffraction  coefficient  and  attenuation  constant  can  be  found  in 
Reference  [47]. 


The  diffraction  thus  far  discussed  is  applied  to  the  open 
curved  surface.  For  a closed  surface,  each  surface  ray  mode 
produced  at  Q-j  encircles  the  surface  an  infinite  number  of  times. 
Trie  length  of  the  surface  ray  path  for  the  nth  encirclement  is 
t+nT  where  T is  the  circumference  of  the  closed  surface.  These 
multiple-encircling  rays  may  be  summed  to  contribute 


rT 


-jkT  - 


am(f)  dt' 


to  the  denominator  of  the  diffracted  field.  It  is  interesting  to 
note  that  there  must  be  another  pair  of  diffraction  points,  Qq  and 
jC,  for  the  closed  surface  as  shown  in  Fig.  4.  Therefore,  the 
field  at  any  point  P in  the  shadow  region  is  the  sum  of  these  two 
a.ffracted  fields  from  Q]-Q2  and  Q3-Q4.  The  total  field  at  any 
point  in  the  illuminated  region  is,  by  the  superposition  principle, 
tne  sum  of  incident,  reflected  and  diffracted  fields.  A detailed 
discussion  of  this  subject  can  be  found  in  Reference  [47], 


3.  DIFFRACTION  BY  AN  INFINITELY  LONG  ELLIPTIC  CYLINDER 

An  important  special  case  of  this  GTD  solution  is  the  one  in 
wnich  the  antenna  is  mounted  directly  on  the  curved  surface.  This 
is  especially  true  for  our  purpose,  since  the  radiation  from  slots 
and  monopoles  mounted  on  smooth  curved  surfaces  is  pertinent  to 
the  design  of  flush-mounted  antennas  for  aircraft  and  spacecraft. 
Recently,  Pathak  and  Kouyoumjian  [48-51]  have  extended  the  GTD 
technique  to  treat  the  radiation  from  apertures  or  slots  in  convex 
perfectly  conducting  surfaces.  This  extension  of  GTD  has  been 
successfully  applied  to  circular  and  elliptic  cylinders,  spheres, 
and  spheroids.  A similar  GTD  analysis  of  the  radiation  from 
monopoles  on  a convex  surface  has,  also,  been  recently  completed 
by  Pathak  and  Luebbers  [52]  with  the  same  degree  of  success. 
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Using  the  GTD  solution,  a launching  coefficient  is  intro- 
duced to  relate  the  antenna  field  to  the  boundary  layer  surface 
waves  which  propagate  around  the  surface  along  geodesic  paths. 
Energy  is  continually  diffracted  by  the  surface  wave  in  the 
tangent  direction  to  the  propagation  path.  This  diffracted 
energy  is  related  to  the  surface  rays  by  a diffraction  coef- 
ficient which  is  dependent  on  the  surface  geometry  at  the  point 
of  diffraction.  The  surface  wave  energy  decays  alone  the 
geodesic  path  in  that  energy  is  continually  diffracted.  This 
decay  is  expressed  by  an  attenuation  coefficient  which  is 
dependent  on  the  surface  geometry  along  the  geodesic  path. 

The  GTD  solutions  for  infinitesimal  slot  and  monopole 
antennas  mounted  on  an  elliptic  cylinder  as  shown  in  Fig.  S 
are  given,  with  torsional  effects  included  [49-92],  by: 

a.  Monopole  Case 

(1)  Illuminated  region 


E = - sin  eme  F(source) 
m m 


(2)  Transition  region 

(a)  Illuminated  side 
n'  x s. 


(31) 


E = 


z x s. 


g*U) 


+ b T (sgn) 


pt  (q') 


■\kpg(Q  ')J 


g*U)  e 


-j  (ks  + C/3) 


f • F (tangent) 


(32) 
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aA 


(b)  Shadow  side 


E = 


^0_| 

1 / 6 

PJL&).  1 

x d1^ 

l g(Q')  J ^ 

n g*(£)  + 


b (sgn)  -1 

J pt(Q') 

e’^^'  • F (tangent) 


L ^ C:pg  ( Q ' ) J\ 


n/3 


g*(0 


(33) 


(3)  Deep  Shadow  Region 


E - I 
j 


n.  E1?  + bjsgn)!-9-^-^ 

J J J . r.  T 


Fj  (tangent) 


(34) 


b.  Slot  Case 

(1 ) 111 uminated  region 


E = [(e1  sin  B - e^  cos  B)  x s]-  F (source) 


(35) 


(2)  Transition  Region 
(a)  Li t side 


E = \ j 

t 

1 n ^ 

sin  a cos  B - 

n‘  x s 

0 

COS  a sin  B 

2 

i 

s 

Iz  X sol 

s 

g*(0 


Pt(Q’ ) 


• so 

C r 

_ 2 " 

9 

1/3 

c 1 

- kpg(Q' ) - 

sin  as  sin  6g*(Of  e 


-j(ksQ  + SJ/3) 


F (tangent) 


(36) 
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(b)  Shadow  side 


,1/6 


-ilar^W-)]  1"  t*i"  K -8)  9*(t)l 


+ b 


M'l  (.  j[  _?  ... 
-Pt(Q' ) ^ |_  k^g(Q* 


e"^  • F (tangent) 


1/3 


^ sin  a$  sin  s§*U) 


(37) 


(3)  Deep  Shadow  Region 


E = \ 


I nj  sin(a$  - 


Fj  (tangent) 


6)  + b.  ^ sin  a sine  E^ 

J J pt(q<)  s j 


(38) 


where 


. | 1 , . -/ ym(n)dft. 

h _ I o r nh  , h 'nr 


E"  = y D"  L"  e 

I d^  '■ n m m 
m-0 


ES  = 


d*0  1 , s -/Ym(Oda 

a—  y L*  e m 
d*  m^0  m m 


sgn 


= t1  • (n1  x z1) 
It*  • (n1  x z 1 ) | 


a + 


1/3 


Y'Yk|  for  illumi"ated  Side 

j q • y 2 Pg 2(i)J  shadow  side 


g*(-),g*(-)  complex  conjugates  of  the  Fock  functions  [48] 
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n' ,b' ,t' ,n,b,t  normal,  binormal,  and  tangent  unit 

vectors  to  the  surface,  at  source  point  Q' 
and  diffraction  point  Q,  respectively. 


F(  • ) 


phase  factor  required  to  refer  the  phase 
of  a ray  to  the  center  of  the  coordinate 
system 

spread  factor  and  equals  unity  for  this 
case  [48] 


p ,p  longitudinal  and  transverse  radii  of 

" T curvature 

z',t'  unit  vectors  in  the  principal  directions 

in  which  the  geodesics  have  no  torsion 

sQ  unit  vector  pointed  from  the  source  to 

the  observation  point 


s the  distance  from  the  source  to  the 

observation  point 


Note  that  the  superscripts  h and  s indicate  the  hard  and  soft 
boundary  conditions,  respectively. 


The  launching  coefficients  are  given  by  [48] 


L 


h 

m 


ire 


JU/12) 


at  the  source  Q' 


L 


s 

m 


at  the  source  Q' 


where  9m  is  defined  in  Table  1.  The  subscript  m refers  to  the 
m*n  mode  of  the  boundary  layer  surface  wave.  Thus,  ym  is  the 
propagation  constant  for  the  mth  mode  surface  wave  such  that 
Ym  = '■•‘m  + jk»  where  um  is  defined  in  Table  I.  The  quantity 
q-,  is  the  root  of  Miller-type  Airy  function  A-j  ( -qm) » that  is, 
Ai(*Pm)*0-  Similarly  qm  is  the  root  of  A|(-qm),  the  derivative 
of  Aj(-qm)  with  respect  to  the  argument  of  the  Airy  function. 
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The  incremental  arc  length  along  the  geodesic  path  is  expressed 
by  di.  The  summation  over  j in  the  shadow  region  indicates  that 
several  terms  can  contribute  in  that  region.  It  is  noted  that  for 
a slot  oriented  at  45°  with  respect  to  the  z-axis,  Eqs.  (36)  to 
(38)  are  not  sufficiently  accurate  for  calculating  the  radiation 
pattern  as  a function  of  <j>  where  e = 46°.  For  this  special  case 
and  for  cases  sufficiently  close  to  it,  an  additional  correction 

term  Ec  should  be  included  as  indicated  in  Reference  [50,51]. 

One  must  first  find  an  efficient  solution  for  the  geodesic 
paths  on  the  elliptic  cylinder  surface  in  order  to  analyze  this 
problem  successfully  using  GTD.  A preferred  coordinate  system 
for  the  elliptic  cylinder  is  illustrated  in  Fig.  6 and  defined  by 

x = d cosh  u cos  v = a^  cos  v 

y = d sinh  u sin  v = sin  v (39) 

z = z 

where  2d  is  the  distance  between  the  foci  of  the  ellipse. 

Note  that  for  u = Uf,  where  Uf  = tanh'^ (bf/af) , the  preceding 
equations  define  an  elliptical  surface.  Thus,  any  point  on 
the  elliptical  surface  is  expressed  by  v,  which  varies  from 
0 to  2n. 

Using  the  calculus  of  variations,  the  z-coordinate  of  the 
geodesic  paths  on  an  elliptical  surface  are  given  by 


2.2  22 
af  sin  v + bf  cos  v dv  + z(v-) 

(40) 

Note  that  Vj  and  Vf  are,  respectively,  the  initial  and  final 
values  of  v along  a given  geodesic  path.  If  one  defines  the 
geodesic  starting  direction  by  the  angle  (ns)  as  shown  in 
Fig.  5,  then  C = -cos  as.  The  advantage  of  this  geodesic 
solution  lies  in  the  fact  that  the  integral  can  be  quickly 
evaluated  using  numerical  techniques.  The  important  parameters 
of  this  problem  are  listed  below: 
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Figure  6.  Diagram  showing  the  elliptic  cylinder  coordinate  system. 
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~~A~ 


(curval inear  coordinates) 


e-,  = 


-a^  sin  vx  + b^  cos  vy 

I 2 7~ 2 2 2 

J af  sin  v + bf  cos 


e2  = 2 

t = sin  ase1  - cos 


(unit  tangent  vector) 


n = 


b ^ cos  vx  + a^  sin  vy 


2.2  22 
a f sin  v + b^  cos 


(unit  normal  vector) 


b = t x n = -cos  a$e1  - sin 


(af2  sin2  v + bf2  cos2  v)3/2 


a£br  sin  a 
f f s 


(unit  bi normal  vector) 


(longitudinal  radius  of  curvature) 


Using  the  preceding  relations,  one  can  employ  ( 31 ) - ( 38)  to 
determine  the  total  radiated  fields  in  the  whole  space  except  for 
two  small  sectors  around  the  cylinder  axis  where  as  is  near  0 or  v , 
since  the  solution  fails  in  these  regions. 


4.  NEAR  FIELD  SCATTERING  BY  A FINITE  BENT  PLATE 

The  near  field  scattering  by  a finite  bent  plate  is  a 
relatively  new  topic  at  higher  frequencies  where  the  plate  is 
large  in  terms  of  the  wavelength.  The  solution  presented  here 
is  a practical  application  of  the  three-dimensional  wedge  dif- 
fraction theory  given  earlier.  The  source  is  defined  by  its 
location  and  far-field  pattern.  The  far-field  pattern  of  the 
source  is  appropriate  in  that  the  plate  is  located  at  least 
2D 2/ > away  from  the  source  where  D is  the  maximum  dimension  of 
the  source.  The  finite  plate  is  simply  specified  by  location 
of  its  n corners.  The  junction  edge  is  defined  by  the  first 
corner  specified  plus  an  additional  corner  (MC)  defined  as  input 
to  the  computer  program.  The  plate  is  initially  flat  (a  = 180°  in 
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Fig.  7).  It  can  then  be  bent  about  the  line  joining  corner  #1 
and  MC  such  that  90°  < a < 270°. 


X SOURCE 
LOCATION 

(*•  ,y$  .**> 


MC+2 


STATIONARY 


PLATE 


Figure  7.  Bent  plate  geometry. 

It  is  known  that  for  a given  scatter  direction  there  is  only 
one  point  along  an  infinitely  long  straight  edge  at  which  the 
diffracted  field  can  emanate  for  a near  zone  source.  Thus,  this 
point  must  be  found  for  each  of  the  n edges  that  describe  the 


finite  plate.  There  are  many  ways  of  finding  this  diffraction 
point,  one  of  which  is  described  here.  Since  it  is  known  that 
S0  = 60‘  (see  Fig.  1),  it  is  obvious  that 


e 


(41) 


where  em,  I,  and  d are,  respectively,  the  mt*1  edge  unit  vector, 
incident  direction  unit  vector,  and  diffraction  direction  unit 
vector.  Since  the  scatter  direction  is  known  (es»$s)»  the 
value  of  em  ■ d = cm  is  easily  computed  at  each  edge.  One  needs 
only  search  along  the  edge  to  find  the  point  where  e • I = c^. 

Once  the  diffraction  point  is  located,  one  must  find  the 
diffracted  field  value  from  the  mth  edge.  The  far  field  pattern 
of  the  source  can  be  written  as 


,-.jks 1 


Es  ( Q »❖ ) = [9  F(9,<tO  + 4>G(0,<f)]  — $1 


= R(e  ,4>) 


(42) 


Vvsere  s'  is  the  range  from  the  source  to  the  field  point.  Using 
tr.e  geometry  illustrated  in  Fig.  7 and  applying  the  results 
presented  earlier,  one  finds  that 


Kll 


-vr 


L o 


lR1 


-j [k(s  * 


y)  -kpp] 


(43) 


wnere 

«i  ■ • v 

. > 

Rl  = R(&i  .4^)  • 4>Q 

2 

k p - ks1  sin  B 
p o 

y = xdp  sin  os  cos  <f$  + ydp  sin  o$  sin  + zdp  cos  e$ 

VB  = VB(kpp/k  ’ 2)  * VB(kpp/k  * *+*‘*  2) 
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The  coordinates  (xdp»  ydp,  zdp)  define  the  point  of  diffraction. 

Once  these  terms  are  determined,  the  total  diffracted  field 
in  ray  form  from  a general  mth  edge  is  given  by 


s’ 


(44) 


wnere  $ = d x B0.  Using  the  superposition  principle,  the  total 
singiy  diffracted  field  in  ray  form  by  the  n edges  of  the  plate 
is  given  using  Eq.  (44)  by 


Rd(e 


s’ 


♦s>  = 


n 

I 

m=l 


Rd(9  . 

nr  s 


♦s) 


(45) 


The  first  step  in  calculating  the  reflected  field  is  to  find 
the  locations  of  the  image  sources,  which  are  uniquely  determined 
once  the  planes  of  the  flat  plates  are  defined  relative  to  the 
source  location.  In  fact,  the  image  is  located  along  lines  which 
are  orthogonal  to  the  different  portions  of  the  plate  and  posi- 
tioned an  equal  distance  on  the  opposite  side  of  the  plate. 

With  the  image  position  known,  one  needs  to  determine  if  the 
reflected  field  contributes  to  the  total  scattered  field  using 
the  geometrical  optics  approach.  If  the  reflected  field  is  a 
contributor,  the  ray  from  the  image  source  in  the  scatter  direc- 
tion ( o 5 , $s)  must  pass  through  the  finite  plate  limits.  Thus, 
cne  must  find  the  location  of  the  intersection  point  of  this  ray 
ana  the  plane  containing  the  flat  plate.  This  can  easily  be 
accomplished  using  vector  analysis.  One  can,  then,  predict 
within  certain  limits  whether  this  intersection  point  falls 
within  the  bounds  of  the  finite  flat  plate. 

If  reflections  do  occur,  the  reflected  field  from  the  image 
source  can  be  written  in  ray  form  as 

i^(9s.  *s)  = [er  Fr(es,  *s)  + $r  Gr(e$,  *s)] 

jk[xi  sin  e$  cos  <f>s  + y^  sin  e$  sin  <j>$  + z.  cos  e ] 

(46) 

where  9 and  <p  are  related  to  the  inage  source  coordinate  system 
with  the  image  location  defined  by  (x^,  y^  , z,).  The  functions 
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[F  ( \ , cs)  and  Gr ( 0 r,  $5)]  are  found  by  employing  the  boundary 
conditions  on  the  perfectly  conducting  flat  plate.  The  total 
scattered  field  from  the  flat  plate  is,  then,  given  by 

Rs(es,  *s)  = Rd(e$,  *s)  + Rr(e$,  *s)  (47) 

The  four  basic  terms  included  in  the  present  solution  are  as 
follows: 

(a)  single  diffraction  of  incident  field  as  shown  in 
Fig.  8(a) 

(b)  single  reflection  of  incident  field  as  shown  in 
Fig.  8(b) 

(c)  double  reflection  of  incident  field  as  shown  in 

Fig.  8(c) 

(d)  single  diffraction  of  reflected  field  as  shown  in 
Fig.  8(d). 

-acn  of  the  terms  has  been  illustrated  in  a two-dimensional  view 
just  for  simplicity  in  illustrating  the  mechanisms-  whereas,  the 
actual  solution  is  for  the  three-dimensional  geometry.  These 
terns  have  all  been  incorporated  in  a general  bent  plate  solu- 
tion. Note  that  only  those  terms  are  included  which  are  not 
shadowed  by  another  portion  of  the  bent  plate. 

The  accuracy  of  this  solution  is  illustrated  by  the  example 
shown  in  Fig.  9 where  it  compares  very  favorably  with  measured 
results  taken  for  a a/2  dipole  illuminating  a flat  plate.  This 
solution  is  also  compared  with  results  obtained  using  a moment 
metnod  patch  technique  [50]  as  illustrated  in  Figs.  10  and  11. 

Note  that  for  these  comparisons  the  plate  dimensions  are  quite 
small  in  terms  of  the  wavelength,  which  accounts  for  the  small 
discrepancies  in  the  patterns  between  the  two  solutions.  This 
structure  will  be  incorporated  into  the  aircraft  model  of 
Section  III  in  order  to  account  for  bent  wings,  moving  flaps, 
vertical  stabilizer,  etc. 
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Figure  8a.  Single  diffraction  of  incident  field. 
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Figure  8b.  Single  reflection  of  Incident  field. 
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Figure  8c.  Double  reflection  of  incident  field. 
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Figure  8d.  "ingle  diffraction  of  reflected  field. 
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Figure  9.  Ef  radiation  pattern  for  a small  dipole  mounted  above 
a square  plate  for  6 = 90°  and  0 - <p  - 360°  at 
f = 10.436  GHz. 
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SECTION  III 

MATHEMATICAL  MODELING  OF  AIRCRAFT 
FOR  PATTERN  COMPUTATIONS 


In  this  section,  the  mathematical  modeling  of  an  antenna 
on  an  aircraft  is  described.  This  mathematical  model  is  intended 
to  be  as  simple  as  possible  so  that  it  can  be  used  to  simulate  a 
wide  variety  of  aircraft  structures,  yet  accurate  enough  that  the 
computed  radiation  patterns  are  comparable  to  measured  results. 

The  model  is  systematical ly  increased  in  complexity  to  resolve 
discrepancies  between  calculated  and  measured  results,  yet  basic 
simplicity  is  retained.  One  of  the  restrictions  in  the  model  is 
that  the  antenna  under  consideration  is  located  near  the  top  or 
bottom  of  the  fuselage  surface. 

This  section  begins  with  a brief  review  of  the  analyses 
using  the  elevation  plane  model  and  the  roll  plane  model.  This 
is  followed  by  solutions  for  the  volumetric  patterns  of  antennas 
on  a three-dimensional  aircraft  structure.  The  approach  used'  to 
model  the  aircraft  by  two  elliptic  cylinders  is  explained  in  detail. 
Tr.e  numerical  technique  to  generate  the  necessary  elliptic  cylinders 
is  also  discussed.  Finally,  the  radiation  patterns  in  principal 
planes  for  various  cases  are  calculated  to  illustrate  the  versa- 
tility of  the  newly  developed  solutions. 


l.  ANALYSIS  OF  ELEVATION  PLANE  MODEL 

In  most  cases,  the  dominant  structural  effect  in  the 
elevation  plane  is  the  profile  of  the  aircraft  for  fuselage 
mounted  antennas.  An  aircraft  fuselage  is  usually  a convex 
body  tnat  cannot  be  completely  described  by  simple  mathematical 
equations.  Hence,  in  practice,  an  aircraft  shape  is  often 
specified  by  points.  Consequently,  a method  called  "section 
matching  method"  has  been  developed  [28]  to  handle  the  problem  of 
an  antenna  mounted  on  a fuselage  surface  of  general  shape. 

m. owever,  this  method  is  restricted  to  principal  plane  pattern 
calculations.  For  off-principal  plane  pattern  calculations 
tne  method  is  hard  to  apply,  since  a complete  three-dimensional 
description  of  the  fuselage  surface  is  difficult  to  obtain. 

Far therm, ore,  the  lack  of  information  on  the  geodesic  paths  and 
tarsion  effect  of  the  curved  surface  create  serious  problems  in 
tne  prediction  of  antenna  patterns.  Consequently,  a simple 
model  which  simulates  a general  fuselage  profile  has  been  adopted. 
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Since  an  aircraft  fuselage  is  usually  long  and  slender,  a 
composite  elliptic  cylinder  appears  to  be  a good  model.  This 
model  consists  of  two  semi-elliptic  cylinders  mounted  back-to- 
back  as  shown  in  Fig.  12.  Note  that  as  discussed  in  the  last 
chapter,  the  model  employed  here  is  not  limited  to  the  determi- 
nation of  the  elevation  plane  pattern  only.  It  can  be  used  to 
determine  off-elevation  plane  radiation  patterns  also. 

Recall  that  the  solutions  presented  in  the  last  section  are 
based  on  elliptic  cylinder  models;  however,  one  of  the  advantages 
of  GTD  is  that  it  can  be  extended  to  new  structures  after  making 
certain  assumptions.  In  this  case,  it  is  assumed  that  there  are  no 
diffractions  from  the  junction  lines  of  the  two  ellipses.  This 
assumption  is  justified  since  these  junctions  are  non-existent  in  the 
actual  aircraft  profile.  Note  that  the  GTD  solution  in  the  illuminated 
region  does  not  depend  on  the  surface  parameters  in  that  it  is 
assumed  the  source  is  mounted  on  an  infinite  ground  plane  tangent 
to  the  surface  at  the  source  point.  On  the  other  hand,  the  transition 
and  deep-shadow  region  solutions  are  modified  due  to  their  dependence 
upon  the  surface  parameters.  This  modification  requires  that  one 
use  bf  = b1  for  rays  traveling  to  the  right  of  the  junction  and 
bf  = b for  rays  traveling  to  the  left  of  the  junction.  The  parameters 
b and  b'  are  illustrated  in  Fig.  12. 
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Since  the  model  employed  is  a composite  elliptic  cylinder, 
the  surface  curvature  at  the  junction  line  is  not  continuous  due 
to  the  different  dimensions  of  the  elliptic  cylinders.  To  ensure 
the  assumption  made  earlier  that  no  diffraction  occurs  from  the 
junction  line,  the  solution  presented  in  paragraph  3 of  Section 
II  needs  to  be  modified  to  correct  for  the  difference  in  curva- 
ture when  the  surface  wave  crosses  the  junction  line.  Since  no 
GTD  solution  for  this  problem  is  available  at  present,  an  alternate 
way  to  account  for  the  continuity  of  the  surface  wave  crossing  the 
junction  line  is  necessary.  It  is  noted  that  due  to  the  large  radius 
of  curvature  at  the  junction  line  in  our  model  the  reflection 
effect  is  so  weak  that  it  can  be  neglected.  A way  to  treat  the 
continuity  of  the  surface  wave  across  the  junction  is  investigated 
next. 


Consider  a source  located  on  a composite  elliptic  cylinder 
as  shown  in  Fig.  13.  The  field  diffracted  by  a curved  surface 
can  be  written,  according  to  Eq.  (34)  for  hard  boundary  con- 
ditions and  with  only  the  first  mode  retained,  as 


E = CL (Q * ) D(Q) 


e 


a + jk)  ds 


(48) 


where  C is  a complex  constant.  This  result  indicates  that  the  field 
diffracted  from  the  curved  surface  is  related  only  to  the  launching 
coefficient  at  the  source  point  and  the  diffraction  coefficient  at 
the  diffraction  point.  Thus,  when  the  diffraction  point  is  on  the 
right  half  elliptic  cylinder,  the  field  is  given  by 

fQ 

(a  + jk)  ds 

E = C L(a,b)  D2(a,b')  e Q‘  (49) 


note  that  a and  b are  the  semiminor  and  semimajor  axes  of  the  left 
half  ellipse  and  a and  b'  are  the  semiminor  and  semimajor  axes  of 
the  right  half  ellipse,  respectively.  Similarly,  the  field 
diffracted  from  point  Q on  the  left  half  ellipse  is  expressed  as 

- I (a  + jk)  ds 

E = C L(a,b)  O^a.b)  e Q’  (50) 
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At  the  junction,  Eqs.  (49)  and  (50)  become 

rJ 


Ej  = C L(a,b)  D2J(a,b*)  e 


(a  + jk)  ds 


and 


Ej  = C L(a,b)  D-|j(a,b) 


'Q' 


(a  + jk)  ds 


respectively,  where  the  subscript  J denotes  the  junction. 

It  is  noted  that  Eqs.  (51)  and  (52)  do  not  agree  with  each 
ther  which  indicates  a discontinui ty  exists  in  the  diffracted 
.ield.  However,  this  is  not  true  in  practice,  since  the  field 
rrjst  be  continuous  as  the  diffraction  point  crosses  the  junction 
Co;  paring  Eqs.  (51)  and  (52),  one  notes  that  if  a modification 
factor 

r 

D i j (a , b ) 

D^TaTb"^ 

is  introduced,  the  field  is  continuous  as  the  diffraction  point 
crosses  the  junction.  Thus,  Eq.  (48)  is  modified  to  give 


0-|j(a,b) 


E • C D^TFFT  L(a'b)  ¥a'b'>  e 


(a  + jk)  ds 


Notice  that  Eq.  (53)  is  valid  only  for  the  surface  wave  which 
crosses  the  junction,  i.e.,  the  diffraction  point  and  source 
point  occur  on  different  elliptic  cylinders. 

From  Table  I,  one  sees  that  the  factor  introduced  in  Eq. 
(53)  can  be  simplified  to  give 


D] j(a,b) 

^jfa.bT 


3gj(a,b) 

p"9j  TTFT 


1/6 


b 

b' 


1/3 


which  depends  only  upon  the  ratio  of  the  two 
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of  the  composite  ellipse.  For  a continuous  elliptic  cylinder, 
this  factor  reduces  to  unity  as  expected.  Note  that  Eq.  (53) 
can  be  generalized  to  describe  the  field  diffracted  from  an 
arbitrary  point  on  the  curved  surface  as 

fQ  , 

(a  + jk)  ds 

E = C L(Q’)  D(Q)  p e Q'  (55) 


where 


r 


if  Q and  Q'  are  on  same  elliptic  cylinder 
ll/6 


^ g j (el  1 ipse  with  source) 
p,g'j(ell  ipse  without  soured 


otherwise. 


(56) 


In  a similar  manner,  this  factor  can  be  determined  through 
the  field  expression  in  the  transition  region  as  presented  in 
Chapter  II  and  is  found  to  be  the  same  as  Eq.  (56). 


To  illustrate  the  newly  modified  solution,  the  elevation 
plane  radiation  patterns  for  antennas  mounted  on  the  KC-135  air- 
craft are  now  analyzed.  Since  the  antennas  of  interest  are  located 
on  top  of  the  fuselage  and  along  the  center  line,  the  most  significant 
effects  on  the  pattern  result  from  the  upper  surface  of  the 
fuselage  profile.  The  structure  used  to  simulate  the  l/25th 
scale  model  consists  of  a 3.75"  by  60"  right  semi -el  1 i pti c cylinder 
and  an  3.75"  by  3"  left  semielliptic  cylinder  as  shown  in  Fig.  12. 
-sing  Eqs.  (31)  to  (38)  in  Section  II  with  the  modifications  just 
described,  the  radiation  patterns  in  the  elevation  plane  for 
antennas  mounted  on  our  model  can  be  easily  determined. 

The  elevation  plane  patterns  for  a \/4  monopole  mounted 
forward  and  over  the  wings  are  illustrated  in  Fig.  14  (a)  and  (b), 
respectively.  The  patterns  for  a circumferential  KA-band  wave- 
guide are  illustrated  in  Figs.  15  (a)  and  (b).  Finally,  the 
results  for  an  axial  KA-band  waveguide  are  shown  in  Fig.  16  (a) 
and  (b).  The  KA-band  waveguide  aperture  fields  are  simulated  in 
our  model  by  an  array  of  infinitesimal  elements  as  shown  in  Fig. 

17. 
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180° 


Figure  14a.  Elevation  plane  pattern  for  a A/4  monopole  mounted 
forward  of  the  wings  on  a KC-135  aircraft.  (No 
radome  and  vertical  stabilizer  included). 


Figure  14b.  Elevation  plane  pattern  for  a a/4  monopole  mounted  above 
the  wings  on  a KC-135  aircraft.  (No  radome  and  vertical 
stabilizer  included). 
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Figure  15a.  Elevation  plane  pattern  for  a ci rcumferentiul  KA-b and 
waveguide  mounted  forward  of  the  wings  on  a KC-135 
aircraft.  (No  radome  an?  vertical  stabilizer  included). 


180° 


Figure  15b.  Elevation  plane  pattern  for  a circumferential  KA-ba.nd 
waveguide  mounted  above  the  wings  on  a KC-135  aircraft. 
(No  radome  and  vertical  stabilizer  included). 
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0 = 0° 


Figure  16a.  Elevation  plane  pattern  for  an  axial  KA-band  wave- 
guide mounted  forward  of  the  wings  on  a KC- 135  air- 
craft. (No  radome  and  vertical  stabilizer  included). 


Figure  16b.  Elevation  plane  pattern  for  an  axial  KA-band  wave- 
guide mounted  above  the  wings  on  a KC- 135  aircraft. 
(No  radome  and  vertical  stabilizer  included). 
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Figure  17.  Approximation  of  an  open-end  waveguide  by  an  array 
of  5 infini tesinal  broad  wall  elements  with  dif- 
ferent weights  and  an  array  of  3 infinitesimal 
narrow  wall  elements  with  uniform  weight. 
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The  comparison  between  calculated  and  measured  results  is 
encouraging,  especially  wnen  the  antennas  are  located  above 
the  wings.  The  discrepancy  displayed  in  the  fore  and  aft  regions, 
particularly  for  the  case  of  antennas  located  forward  of  the  wing, 
is  due  to  the  effects  of  the  cockpit  and  nose  section  and  the 
vertical  stabilizer  which  are  neglected  in  the  present  theoretical 
model . 


In  the  preceding  analysis,  the  effect  of  the  aircraft  in  the 
nose  of  the  radome  was  neglected.  In  order  to  simulate  a fuselage 
with  a radome,  a truncated  composite  elliptic  cylinder  model  as 
snown  in  Fig.  13  is  adopted.  This  is  based  on  experimental  data 
that  the  radome  effect  appears  to  be  a diffraction  from  the  bulk 
head  as  if  the  radome  were  not  present. 

Before  proceding  to  solve  this  problem,  by  careful  consider- 
ation it  can  be  converted  into  a simpler  problem.  With  the  fact 
tnat  a slot  or  monopole  radiating  in  the  presence  of  a composite 
elliptic  cylinder  can  be  considered  as  an  antenna  itself,  the 
present  problem  can  be  approximated  by  the  radiation  problem  of 
an  equivalent  antenna  mounted  on  a wedge  type  structure  as  shown 
in  Fig.  19. where  the  antenna  is  mounted  on  the  bottom  of  the 
fusel >ge.  This  equivalent  antenna  radiates  the  same  pattern  as 
toet  of  a slot  or  monopole  mounted  on  our  previous  model  without 
toe  radome  considered.  The  wedges  are  formed  and  defined  by  the 
tangent  planes  at  the  discontinuities  0]  and  Q^.  Thus,  the 
radiation  problem  of  a fuselage  model  with  a radome  included 
is  transformed  to  a wedge  diffraction  problem,  and  can  be  solved 
using  standard  GTD  techniques. 
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Figure  18.  Truncated  composite  ellipse  as  simulation  of  fuselage 
with  radome  included. 


Figure  19.  Total  field  distribution. 


Using  Eq.  (57),  the  elevation  plane  patterns  for  the  space 
shuttle  as  shown  in  Fig.  20  with  the  radome  included  are  calcu- 
lated. Figures  21  to  23  illustrate  the  radiation  patterns  for  both 
vertical  (monopole  and  circumferential  waveguide)  and  horizontal 
(axial  waveguide)  polarization  sources  mounted  on  the  fuselage  of 
the  1 / 35th  scale  model  of  the  space  shuttle  with  the  radome  considered. 
The  structure  used  to  simulate  the  1 / 35th  scale  model  consists 
of  a 2.55"  by  12.5"  left  semi-ellipse  and  a 2.55"  by  60.0"  right 
semi-ellipse.  The  size  of  the  radome  is  6.5"  and  the  source  is 
mounted  2"  away  from  the  radome.  Figure  24  shows  the  radiation 
pattern  of  the  space  shuttle  with  radome  size  being  2"  and 
antenna  mounted  1 fZ"  away  from  radome.  Mote  that  the  waveguide 
used  here  is  the  same  KA-band  waveguide  used  previously.  The 
experimental  results,  which  were  taken  at  NASA  (Hampton,  Va.)  at 
a frequency  of  35  GHz,  are  also  presented.  The  agreement  between 
calculated  and  measured  results  illustrates  the  applicability  of 
the  numerical  analysis  technique  being  developed  in  this  study. 

Mote  also  that  although  the  elevation  plane  model  can  be  used 
to  compute  the  elevation  plane  pattern  as  well  as  off-elevation 
plane  patterns,  it  is  found  that  for  our  fuselage  model,  this 
model  can  only  be  used  to  compute  the  off-elevation  plane  patterns 
accurately  up  to  approximately  -30°  from  the  elevation  plane. 

2.  ANALYSIS  OF  ROLL  PLANE  MODEL 

The  basic  aircraft  to  be  analyzed  in  this  section  is  composed 
of  flat  plates  attached  to  an  infinitely  long  elliptic  cylinder. 
Cricinally,  a circular  cylinder  was  used  to  represent  the  fuselage 
[27j,  and  the  modal  solutions  were  employed  to  determine  the 
field.  However,  an  aircraft  cross  section  is  not  circular  in 
general.  To  be  able  to  predict  the  radiation  pattern  more 
accurately,  one  has  to  consider  a better  model  for  the  fuselage 
approximation. 

Since  the  roll  plane  cuts  orthogonally  across  the  fuselage, 
cne  should  expect  the  fuselage  cross  section  to  have  a strong 
effect  on  the  roll  plane  pattern.  On  the  other  hand,  an  aircraft 
fuselage  is  normally  long  and  slender,  such  that  its  finite  length 
e*fects  are  generally  secondary.  Consequently,  the  infinitely  long 
elliptic  cylinder  representation  of  the  fuselage  for  roll  plane 
calculations  appears  to  be  a reasonable  approximation  in  most 
cases.  Since  tne  antenna  can  be  arbitrarily  positioned  on  the 
fjselage  with  respect  to  the  wings,  one  must  consider  the  width 
of  tne  wing  as  well  as  its  length  in  order  to  obtain  a practical 
analytic  model.  In  order  to  accomplish  this,  the  near  field 
oent  plate  scattering  solution  is  adapted  to  this  new  model 
such  as  illustrated  in  Fig.  25.  Note  that  the  wings  are 
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Figure  20a.  Scale  model  (1/35)  of  Space  Shuttle  in  Anechoic 
Chamber  at  NASA,  (Hampton,  Va.). 
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Figure  20b.  Theoretical  model  of  1 / 3 5 th  scale  model  of  space 
shuttle. 
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VERTICAL  POL.  SPACE  SHUTTLE 
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Figure  21.  Elevation  pattern  of  monopole  on  a l/35th  scale  model 
of  space  shuttle  (with  radome  included). 
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Figure  22.  Elevation  pattern  of  ci rcumferential  slot  on  a 
l/35th  scale  model  of  space  shuttle  (with  radome 
included). 
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Figure  23.  Elevation  pattern  of  axial  slot  on  a 1 /35th  scale 
model  of  space  shuttle  (with  radome  included). 
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Figure  24.  Elevation  pattern  of  axial  slot  on  a l/35th  scale 
model  of  space  shuttle  (with  radome  being  2"). 
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Figure  25a.  Reflection  problem  in  x-y  plane. 


Figure  25b.  Diffraction  problem  in  y-z  plane. 


58 


assumed  flat  here  and  each  wing  can  be  located  arbitrarily 
with  any  number  of  edges.  The  wings,  also,  can  be  mounted  up 
or  down  from  the  central  location  provided  that  the  wings  are 
horizontal . 


The  model,  now,  consists  of  an  infinitely  long  elliptic 
cylinder  fuselage  to  which  finite  flat  wings  are  attached.  The 
various  configurations  analyzed  are  shown  in  Fig.  26  looking 
from  the  front  of  the  aircraft  with  the  antenna  mounted  in 
each  case  above  the  wings  for  the  models  illustrated.  Using 
tnese  models,  one  should  be  able  to  analyze  a wide  variety  of 
aircraft  shapes.  This  is  verified  by  a comparison  of  results 
taken  on  actual  aircraft  scale  models  and  presented  later. 

Let  us  first  find  the  effective  source  location  for  the 
reflected  field.  Recall  that  in  the  flat  plate  result,  the 
source  was  imaged  and  the  reflected  field  added  to  the  total 
solution  provided  the  image  ray  passed  through  the  finite 
flat  plate  (wing)  limits.  One  must  initially  determine  the 
effective  source  position  and  then  the  reflected  field.  With 
the  source  mounted  on  an  infinitely  long  elliptic  cylinder, 
the  surface  rays  from  the  source  propagate  around  the  cylinder 
along  geodesic  paths,  from  which  energy  is  continually  diffracted 
tangentially.  Now  let  us  assume  that  the  source  does  not  illumi- 
nate the  right  wing  directly  (as  illustrated  in  Fig.  25(a))  and 
proceed  to  determine  the  unique  geodesic  path  that  diffracts 
energy  from  a known  tangent  point  which  is  then  reflected  off 
tne  wing  in  the  desired  radiation  (or  scatter)  direction. 


The  effective  source  position  for  reflections  from  the  right 


wing  in  terms  of  the  radiation  direction  (or,  |S)  is  given  by 


= df  cos  v, 
e f e 


y - bf  sin  v 
Je  f e 


(58) 


z = cot  0 
e s 


•1 


so 


2.2  2 2 

a,  sin  v + br  cos  v dv  + z 
f f so 


where  ve  = tan  ( bf /a f cot  <1^).  These  coordinates  car , then, 
be  used  in  the  flat  plate  problem  as  the  effective  source  loca- 
tion. Note  that  as  the  desired  radiation  direction  is  varied, 
the  effective  source  location  changes.  In  addition,  if  the 
source  directly  illuminates  the  wing  for  a given  reflection 
term,  then  the  effective  source  location  is  simply  the  actual 
source  location  (uf,  vs_,  z.lCf) . A result  similar  to  Eq.  (58) 
can  be  found  for  the  reflections  from  the  left  wing.  Finally, 
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Figure  26.  Fuselage  and  wing  geometry  for  theoretical  aircraft 
model  looking  from  the  front.  The  antenna  is  always 
mounted  on  top  of  the  models. 
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the  actual  source  field  value  used  to  compute  the  reflected 
term  is  determined  from  the  GTD  solutions  in  Chapter  II. 


Using  a similar  technique,  the  effective  source  locations 
for  the  diffracted  field  components  may  be  found.  The  bent 
plate  solution  uses  a search  technique  to  find  the  diffraction 
point  by  computing  the  diffraction  angles  at  selected  test 
points  along  a given  edge.  Once  a test  point  (xd,  yd,  zd)  is 
specified  along  the  edge,  one  can  find  the  effective  source 
location  (xe,  ye.  ze)  using  the  geometry  illustrated  in  Fig. 
25(b).  Again  it  is  assumed  that  the  source  does  not  directly 
illuminate  the  test  point.  One  finds  that  the  effective  source 
is  given  by 


x 


e 


2.2 

afbfxd 


+ a 


2 rr 

fydV  af> 


dfy 


.2  2 

bfxd 


,2  2 
Vd 


2,2 

- afbf 
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afbfyd 
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2 f 2 2 A .2  2 2.2 

fxd  1 a*yj  + b-cXj  ‘ a-cb 


f_yd  ■ "fAd 
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afyd  + bfxd 


fuf 


ze  = 


b.x  z.l  + a , z (y , - y ) I ' 
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wnere 


rv  - 


rv* 


V. 


I 2 _ 2 2 2 

J df  sin  v + b^  cos  v dv. 


so 


I' 

V 


'yjb. 


f~2  2 2 2 .i/J'P,uf 

Jdf  sin  vp  + bf  cos  vn,  and  v„  = tan  ' 


*e/a(v 


Given  the  effective  source  location  for  the  chosen  test 
point,  the  search  technique  is  applied  to  find  the  actual  dif- 
fraction point  along  a given  edge.  Note  that  once  the  actual 
diffraction  point  is  determined,  the  effective  source  of  the 
diffracted  field  is  specified  by  Eq.  (59),  and  the  source  field 
value  is,  again,  computed  using  the  GTD  solutions. 
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The  total  field  is  found  by  summing  the  directly  radiated 
field  with  the  scattered  field  using  the  superposition  principle. 
The  results  for  ,j  a/4  inonopole  on  the  fuselage  of  a 1 /25th  scale 
model  of  a KC-135  aircraft,  as  shown  in  Fig.  12,  forward  and  over 
the  wings  are  shown  in  Figs.  27  (a)  and  (b),  respectively.  The 
results  for  a KA-band  circumferential  waveguide  forward  and  over 
the  wings  are  shown  in  Fig.  28.  The  results  for  a KA-band  axial 
waveguide  forward  and  above  the  wings  are  shown  in  Fig.  29.  The 
fuselage,  in  this  case,  is  approximated  by  a 3"  x 3"  infinitely 
long  elliptic  cylinder.  The  waveguide  antenna  is  modeled  as  in 
the  previous  section,  and  the  agreement  in  each  case  is  very 
encouraging.  Again,  it  should  be  stressed  that  the  roll  plane 
model  solution  not  only  can  predict  the  radiation  pattern 
accurately  in  the  roll  plane,  it  also  can  be  extended  to  cover 
almost  tne  complete  volumetric  pattern  as  shown  in  Ref.  [30]. 


3.  THREE-DIMENSIONAL  MODEL  APPROXIMATION  OF  AN  AIRCRAFT 

’’he  analytical  solutions  described  in  the  previous  sections 
provide  a useful,  efficient,  and  economical  way  for  toe  evalu- 
atis^,  Icca'ion , and  design  of  fuselage  mounted  antennas  based 
or  treir  pattern  performance  in  principal  planes.  However,  if 
~cce'-n  systems  are  to  function  properly,  the  antenna  pattern 
".st  meet  certain  specifications.  These  specifications  are 
Uiuaily  given  in  terms  of  a coverage  diagram  for  a particular 
sector  in  space  which  is  to  be  met  with  antenna  mounted  on  the 
air  raft  not  on  a finite  ground  plane  where  the  original  antenna 
det  gr.  war  made.  Thus,  the  desire  for  an  accurate  solution  for  the 
complete  pattern  performance  of  antennas  mounted  on  a complex  air- 
crart  structure  for  given  applications  requires  a more  thorough  study 
of  ways  to  handle  the  volumetric  pattern. 

If  this  problem  is  attacked  directly  by  analyzing  rays 
on  corple/^three-dimensional  surfaces  as  done  previously  in 
Reference  [2],  tne  resulting  numerical  solution  would  be  very 
complex,  time-consuming,  and  uneconomical.  Nevertheless,  if 
certain  assumptions  can  be  made,  the  approach  undertaken  in 
tne  previous  sections  can  be  used  to  overcome  these  difficulties 
and  simplify  the  problem  a great  deal. 

First,  it  has  been  shown  by  numerous  scale  model  measure- 
ments that  the  roll  plane  model  can  be  extended  to  almost 
c.ver  the  complete  volumetric  pattern  except  for  two  conical 
s s'.ors  (fore  and  aft)  [30].  The  limitations  of  the  roll  plane 
model  result  are  due  to  the  finite  length  fuselage.  Yet,  the 
finite  length  fuselage  has  been  solved,  previously,  in  the  elevation 
plane  model  analyses.  Furthermore,  based  on  the  previous 
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'igure  27a.  Roll  plane  pattern  (E.)  for  a 1/25  scale  model  of  a 
KC-135  with  a //4  mon6pole  on  the  fuselage  forward 
of  the  wings  at  freq.  = 34.92  GHz  (model  frequency). 
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Fiqure  27b.  Roll  plane  pattern  (Ej  for  a a/4  monopole  above 
the  wings.  * 
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riqure  28a.  Roll  plane  pattern  (E,  ) for  a KA-band  circumferential 
waveguide  forward  of  {he  wings. 


180° 

Figure  28b.  Roll  plane  pattern  (L,  ) for  a KA-band  circumferential 
waveguide  above  the  wings. 


64 


three-dimensional  studies  of  geodesic  rays  which  contribute  to  the 
pattern  of  an  antenna  on  various  prolate  spheroids,  one  is 
able  to  combine  the  analyses  of  these  two  principal  plane  models 
to  give  the  complete  pattern.  In  this  study,  computer  simulation 
models  were  considered  that  would  resemble  a wide  variety  of 
aircraft  shapes  and  yet  could,  also,  be  analyzed  with  reason- 
able accuracy  and  economy.  In  this  case,  it  is  quite  obvious 
that  the  three-dimensional  nature  of  the  fuselage  must  be 
modeled  if  one  is  to  adequately  determine  volumetric  patterns. 

This  resulted  in  the  development  of  a general  surface  of 
revolution  model  for  the  fuselage  as  presented  in  Reference 
[2].  Through  an  extensive  study  of  geodesic  paths  on  a surface 
of  revolution,  the  number  of  dominant  rays  that  contributed  to 
tne  radiation  pattern  were  shown  to  be  finite  except  for  more 
spnerical  shapes.  Furthermore,  the  computer  result  showed  that, 
for  a prolate  spheroid,  the  dominant  rays  needed  to  be  con- 
S'dered  would  not  exceed  four  rays;  in  most  cases,  it  is  even 
less  than  that.  These  four  rays  are  illustrated  in  Fig.  30 
in  woich  two  rays  are  propagating  along  the  cross-section  of 
prolate  spheroid;  the  other  two  are  propagating  along  the  profile. 


Fioure  30,  The  four  dominant  GTD  terms  that  radiate  at 
e = 90°,  « = 145°). 
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To  demonstrate  the  significance  of  the  four-ray  contribution, 
the  elevation  plane  pattern  of  an  axial  slot  mounted  on  a prolate 
spseroid  was  calculated  using  a two-dimensional  (two  rays)  and 
three-dimensional  (four  rays)  solution  as  shown  in  Fig.  31. 
Experimental  results  are  also  shown  to  verify  the  calculated 
three-dimensional  result.  It  is  immediately  obvious  that  the 
tack  lobe  region  is  not  calculated  with  sufficient  accuracy 
using  the  two-dimensional  result.  However,  the  three-dimensional 
solution  is  in  good  agreement  with  the  measured  pattern.  This 
leads  to  a new  approach  to  handle  the  volumetric  patterns  for 
fuselage  mounted  airborne  antennas  in  a simplified  and  economical 
manner. 


As  determined  previously,  the  roll  plane  model  (infinite 
elliptic  cylinder  with  flat  plate  wings)  can  be  extended  to  cover 
most  of  the  volumetric  pattern  except  for  the  fore  and  aft  sectors 
(about  15°  conical  sector  as  shown  in  Fig.  32).  To  cover  these 
sectors,  one  must  incorporate  the  elevation  plane  model  analysis 
with  some  modifications  such  that  the  effects  of  wings  and  stabilizers 
are  also  included.  In  doing  so,  two  different  elliptic  cylinders  are 
required  to  analyze  the  volumetric  patterns;  one  being  the  cross- 
section  (roll)  cylinder  and  the  other  being  the  profile  (elevation) 
cylinder.  This  requirement  is  necessary  since  our  solutions  are 
based  on  an  infinitely  long  elliptic  cylinder,  in  which  the  z-axis 
coincides  with  the  axis  of  the  cylinder. 


As  discussed  in  the  previous  chanter,  the  most  significant 
effects  on  the  radiation  pattern  result  from  the  surface 
geometry  nearest  to  the  antenna.  The  curvature  of  the  surface 
in  the  vicinity  of  the  antenna  location  plays  a dominant  role 
in  predicting  the  radiation  pattern.  Thus,  the  elliptic 
cylinders  used  to  represent  the  fuselage  profile  and  cross- 
section  need  to  model  the  aircraft  structure  as  accurately 
as  possible  near  the  antenna  location.  Once  these  two 
elliptic  cylinders  are  obtained,  one  is  able  to  proceed  to 
solve  for  the  complete  volumetric  pattern. 


To  determine  the  elliptic  cylinders  necessary  to  simulate 
the  fuselage  profile  and  cross-section,  one  has  to  find  some 
way  to  approximate  the  curved  surface  near  the  antenna  location. 
Ere  way  to  achieve  this  goal  is  to  find  an  analytical  function 


which  can  best  approximate  the  known  curved  surface 
(either  in  elevation  or  roll  plane)  in  a least  mean 
sense.  This  leads  to  the  development  of  a best  fit 
to  generate  the  necessary  ellipses  for  the  fuselage 
through  the  use  of  a digital  computer. 


profile 
square  error 
routine 
model 
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Figure  31a.  Elevation  plane  pattern  of  an  axial  slot  mounted 
on  a 4a  x 2a.  prolate  spheroid  with  the  two- 
dimensional  theoretical  solution  presented. 
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Figure  31b.  Elevation  plane  pattern  of  an  axial  slot  mounted  on 
a 4a  x 2A  prolate  spheroid  with  the  three-dimensional 
theoretical  solution  presented. 
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The  theory  behind  the  best  fit  routine  is  that  a function 
is  to  be  found  to  best  approximate  a set  of  points  in  a least 
mean  square  error  sense.  For  the  aircraft  model,  a best  fit 
ellipse  is  desired.  The  mathematical  expression  for  an  ellipse 
is 


(60) 


for  its  origin  located  at  (x  = 0,  y = 0).  The  parameters  a and 
b are  semi  major  and  minor  axes  of  an  ellipse,  respectively.  To 
simplify  the  mathematical  expression,  Eq.  (60)  can  be  written  as 


AX  + BY  = 1 (61) 

2 2 1 1 

where  X = x,Y  = y,A  = — and  B = —*  . Let  (x-,  y.), 

a b 

i = 1,  •••  n,  be  n points  on  which  a best  fit  ellipse  is  to  be 
generated.  Substituting  these  points  into  Eq.  (61),  one  obtains 
a set  of  n linear  equations. 

AX]  + BY1  = 1 

ax2  + by2  = 1 


AX  + BY  =1 
n n 

In  matrix  form,  these  n equations  become 


ZC  = I 


wnere 


(62) 
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By  multiplying  both  sides  of  Eq.  (62)  by  Z,  the  transposed 
of  Z,  one  obtains  a simple  2x2  matrix  as  given  by 


or 


% ^ 

Z ZC  = Z I 
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(63) 


which  can  be  simply  solved  for  A and  B in  the  least  mean 
sqjare  sense  [56].  Thus,  the  ellipse  parameters  are  defined. 

Before  the  necessary  ellipse  or  composite  ellipse  can  be 
determined  through  the  numerical  process,  the  data  points  that 
described  the  surface  of  the  profile  or  cross  section  of  the 
a tual  aircraft  fuselage  must  bo  generated.  To  do  so,  a reference 
Cartesian  coordinates  system  is  needed  on  a scale  model  drawing 
of  the  aircraft.  These  coordinates  can  be  best  located  by  align- 
ing ore  of  the  axes  with  the  center  line  of  aircraft  fuselage 
with  the  origin  being  arbitrarily  chosen  according  to  convenience. 
After  the  coordinate  system  is  fixed,  the  positions  of  data  points 
on  the  profile  or  cross-section  can  be  measured  from  the  scale 
model  relative  to  the  reference  origin.  The  data  points  are 
tc.ren  in  such  a way  that  more  points  are  needed  around  the 
antenna  location  and  less  points  away  from  the  source.  This  is 
dje  to  the  fact  that  the  surface  profile  is  dominant  near  the 
antenna  location,  as  described  earlier.  Figure  33  illustrates 
the  way  data  points  are  taken  from  a fuselage  profile.  By  feeding 
these  data  points  into  the  best  fit  routine  and  adjusting  the 
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Figure  33.  Illustration  of  data  points  taken  from  the  scale 
model  aircraft  for  the  determination  of  best  fit 
elliptic  cylinder  using  a digital  computer. 
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origin  of  the  coordinates  in  the  routine,  an  ellipse  is  found  to 
best  fit  these  data  points.  In  the  same  manner,  the  ellipse  to 
approximate  the  cross-section  of  aircraft  fuselage  at  the  antenna 
location  can  be  obtained.  Thus,  the  elliptic  cylinders  necessary 
to  simulate  the  fuselage  profile  and  cross-section  are  determined. 

Now  let  us  consider  how  these  two  finite  elliptic  cylinders 
can  be  used  to  describe  the  volumetric  pattern  for  fuselage  mounted 
airborne  antennas.  First,  a reference  coordinate  system  is  needed 
so  tnat  the  solutions  based  on  our  two-dimensional  analysis  in  two 
different  coordinate  systems  can  be  incorporated  into  a complete 
solution.  Since  an  aircraft  fuselage  is  usually  long  and  slender, 
its  finite  length  effects  are  limited  to  the  small  sectors  off  the 
nose  and  tail.  Consequently,  the  two  rays  which  play  the  most 
significant  role  normally  come  from  the  cross-sectional  elliptic 
Cylinder.  Thus,  it  seems  more  natural  to  have  the  reference 
cocrr,  sate  system  correspond  to  the  roll  plane  model  coordinate 
systt  m,  i.e.,  the  Zrc.f  axis  is  pointed  aft  and  the  XrPf  axis  is 
pointed  vertically  upward  as  shown  in  Fig.  34.  Notice  that  the 
origin  of  the  composite  ellipse  in  the  elevation  plane  is  chosen 
as  the  origin  of  tne  reference  coordinate  system  for  convenience. 

Tecall  that  the  roll  plane  model  is  not  valid  in  the  two 

153  conical  sectors  in  the  fore  and  aft  directions.  In  order  to 

overcome  this  handicap,  a belt  region  around  the  antenna  source 
is  cnosen  such  that  in  this  region  the  solution  is  obtained  using 
elevation  plane  model  analysis.  Everywhere  outside  the  belt  region, 
tne  roll  plane  model  analysis  is  used  to  analyze  the  radiation  pat- 
tern as  done  previously.  Figure  35  shows  the  regions  in  which  the 
elevation  and  roll  plane  cylinder  solutions  are  used.  The  angle 
n is  chosen  such  tnat  the  roll  and  elevation  cylinder  solutions 
b ?-d  smoothly  together.  In  fact,  i is  a function  of  the  size  of 
e’evation  ar.d  roll  plane  cylinders  used  to  simulate  the  aircraft 
Tuselage.  for  most  cases  in  our  model,  the  angle  i is  set  at 

'i-j  . This  23°  belt  has  been  tested  and  found  to  be  satisfactory, 

based  on  comparisons  with  measured  results,  as  will  be  shown  later. 

As  seen  in  Fig.  34,  the  elevation  plane  model  coordinate  sys- 
tem is  such  that  the  ZeiGV  axis  points  in  the  Yref  direction.  In 
ctr.'.r  words,  the  elevation  plane  model  coordinate  system  is  simply 
cotained  by  rotating  the  XYZref  coordinate  system  by  a 90°  angle 
about  the  XrGf  axis.  Notice  that  the  antenna  here  is  assumed  to  be 
rGw.rtej  on  the  center  line  of  the  fuselage  for  simplicity,  but  this 
is  rot  a basic  limitation.  In  order  to  determine  the  electromagnetic 
fields  in  a given  radiation  direction  ( o0  ref,  ref)  within  the 
celt  region,  the  solution  obtained  through  the  elevation  plane 
rcdel  analysis  is  expressed  in  terms  of  the  reference  coordinate 
system  such  that  it  can  be  superimposed  with  the  roll  plane 
cylinder  solution.  This  requires  a coordinate  transformation 
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Figure  34a.  Composite  aircraft  model  and  coordinate  systems  illustration 
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* Figure  35.  Illustration  of  belt  region  in  which 

elevation  plane  model  analysis  is 
employed. 
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which  relates  the  elevation  plane  model  coordinates  to  the  reference 
coordinate  system.  The  transformation  of  coordinate  systems  is 
examined  in  Appendix  A. 

Notice  that  the  coordinate  transformations  described  in 
Appendix  A are  very  important  for  the  analysis.  This  is  especially 
true  when  the  antenna  under  consideration  is  mounted  off  the 
center  line  of  the  fuselage.  Using  these  transformations  not  only 
can  the  electromagnetic  field  in  the  desired  radiation  direction  in 
the  reference  coordinate  system  be  determined  through  the  two 
principal  plane  cylinder  solutions,  but  the  polarization  problem 
associated  with  the  two  different  solutions  can  also  be  easily 
resolved  without  any  difficulty.  In  addition,  these  transformations 
are  the  foundation  on  which  the  volumetric  pattern  analysis  is  based. 
A detailed  discussion  on  an  off-center  line  mounted  antenna  is 
given  in  Appendix  B. 

To  determine  the  radiation  angle  (ore,  4>re ) in  the  elevation 
plane  model  coordinate  system  which  corresponds  to  the  desired  radia- 
tion direction  (9y-o,  <£ro)  in  the  reference  coordinate  system, 

Eqs.  (80)  and  (83)  are  employed.  The  subscripts  e and  o refer  to 
the  elevation  plane  model  and  reference  coordinate  systems, 
respectively.  By  definition,  the  unit  radial  vector  for  a given 
direction  (0,  $)  is  given  by 


R = sine  cos<j>  x + sin  0 sin<j>  y + cose  z (64) 


Hence,  the  spherical  coordinates  (0  , <(>  ) of  the  radiation 

vector  are  determined  by  re  re 


-1  y 


re 


♦«,  = tan  — 


re 


re 


re 


tan 


-1 


( X*"  + )^/^ 

're  re' 

Zre 


(65) 


with  the  components  X . and  Z found  through 

re  re  re 
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- 

For  the  case  where  <J>  = 0°  and  0 = 90°  as  seen  in  Fig.  36, 

Eq.  (66)  reduces  to 
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Finally,  the  radiated  field  in  the  desired  radiation  direc- 
tion can  be  determined  from  the  field  found  in  the  elevation  plane 
model  analysis  with  an  appropriate  polarization  modification.  Thus, 
the  solution  obtained  can  be  incorporated  into  the  roll  plane 
cylinder  solution  to  form  a complete  solution  for  the  volumetric 
pattern. 


Let  us  now  consider  the  polarization  effect  due  to  the 
coordinate  transformation,  which  is  vital  in  achieving  a com- 
plete solution.  Recall  that,  in  the  far  zone,  the  electric  field 
can  be  written  as 


E = E0  9 + $ (68) 


or 

E = E-  6 1 + E;  (69) 
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Figure  36.  Illustration  of  the  spherical 
coordinate  system. 


where  the  prime  system  denotes  the  new  coordinate  system  as  seen 
in  Fig.  36.  Since  the  components  Eq  and  El  are  found  from  the 
elevation  plane  solutions,  the  corresponding  E0  and  E*  in  the 
desired  coordinate  system  are  yet  to  be  determined.  If^one 
takes  the  dot  product  of  Eq.  (69)  with  the  unit  vector  0,  the 
e component  of  the  electric  field  results,  i.e. 

E0  = Eg  (01  • e)  + E^  U1  • 0)  (70) 

Similarly,  the  41  component  can  be  found  as 

E^  = Eq  (e*  • t)  + EJ  U1  • i)  (71) 
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Thus,  the  polarization  of  the  radiated  field  in  terms  of  the 
reference  coordinate  system  can  be  properly  determined. 

Before  the  roll  and  elevation  plane  model  solutions  can  be 
combined  together,  one  must  be  concerned  with  the  alignment 
of  tne  surface  normal  at  the  antenna  location  in  the  two 
analyses.  For  an  antenna  mounted  centrally  on  the  fuselage,  the 
surface  normal  points  vertically  unwards  in  the  roll  plane  model 
analysis;  whereas,  in  the  elevation  plane  model  it  points  in  a 
direction  normal  to  the  profile  surface  at  the  antenna  location. 

If  tne  normals  do  not  align  properly,  the  two  solutions  would  not 
match  to  give  a complete  and  continuous  solution.  Consequently, 
tne  roll  plane  elliptic  cylinder  needs  a tilt  such  that  the  normal 
direction  is  common  to  both  solutions.  It  is  noted  that  a radial 
monopole  is  mounted  parallel  to  the  surface  normal  at  the  antenna 
location  on  the  aircraft  fuselage. 


”0  accomplish  this  alignment,  the  Xroll  axis  of  the  roll 
plane  model  coordinate  system  is  chosen  to  line  up  parallel  with 
: .normal  to  the  surface  at  the  antenna  location  on  the  elevation 
cciposite  elliptic  cylinder  as  shown  in  Fig.  34.  The  tilt  angle 
v.,,,  is  found  to  be 

LI  I t 


0 


tilt 


r 

tan 


ae  Zs  ref^ 
be  Xs  ref 


tan 


-II 


' ae  Zs  ref 
^ Xs  ref 


if  Zs  ref  is  ne9dtive 


if  Zs  ref  is  positive 


(72) 


w'ere  Xs  ref  and  Zc  ref  are  the  coordinates  of  the  antenna  source 
relative  to  the  reference  coordinate  system.  The  quantities  ae. 

Dp , and  bg  are  the  necessary  parameters  for  the  composite  elliptic 
cylir.cer  model  for  the  elevation  plane  model  analysis.  Notice  that 
t^e  origin  of  the  roll  plane  cylinders  is  located  a distance  ar 
erectly  under  the  antenna.  Again,  the  quantity  ar  represents  an 
a^.is  of  the  roll  plane  cylinder. 

It  is  also  noted  that,  the  parameter  ar  is  not  equal  to 
Xs  ref  which  is  tne  height  of  the  antenna  on  the  elevation 
composite  cylinder  model.  This  is  due  to  the  fact  that  the 
two  elliptic  cylinders  used  to  approximate  the  fuselage  profile 
and  cross  section  do  not  match  to  give  the  true  representation 
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of  a three-dimensional  aircraft  fuselage.  Hence,  to  obtain  the 
necessary  roll  plane  model  solution,  a coordinate  transformation 
similar  to  that  of  the  elevation  plane  model  analysis  described 
above  is  needed.  The  required  rotation  angle  (o0,  $0)  of  the  roll 
plane  model  coordinate  system  is  found  to  be 


9 _ - 9 , . , . 
0 tilt 


*0  = 0 


(73) 


for  a centrally  mounted  antenna.  Since  the  angle  is 

rotated  about  the  y'-axis,  Eqs.  (80)  and  (85)  shall  be  used. 

One  last  note  for  the  completion  of  our  solution  is  the 
phase  reference  problem.  For  both  elevation  and  roll  plane  model 
analyses,  the  origin  of  the  individual  coordinate  system  is 
chosen  to  be  the  phase  reference  point.  Hence,  the  result 
obtained  through  these  two  principal  plane  model  analyses  needs  to 
be  modified  so  that  a common  phase  center  can  be  utilized.  For 
convenience,  the  origin  of  the  reference  coordinate  system  is 
selected  as  the  phase  reference  point  for  both  models.  In  this 
way,  the  solutions  blend  smoothly  together  to  form  a complete 
result  for  the  three-dimensional  geometry  under  cons ideration. 

For  each  of  the  two  solutions  just  described,  only  two 
rays  are  considered.  This  is  due  to  the  fact  that  there  are 
only  two  dominant  rays  which  propagate  around  an  electrically 
large  elliptic  cylinder.  Recall  that  four  dominant  rays 
contribute  to  the  radiation  pattern  in  the  shadow  region  for  an 
antenna  mounted  on  a prolate  spheriod.  To  account  for  these 
four  rays  necessary  for  the  radiation  pattern  in  the  shadow 
region,  both  elevation  and  roll  plane  model  results  are  needed 
in  the  pattern  calculation. 

In  addition,  the  three-dimensional  effects  of  wings  and 
stabilizers  are  considered  by  including  the  rays  reflected  and 
diffracted  from  these  scatterers.  The  contribution  of  wings 
and  horizontal  stabilizers  are  handled  in  the  same  way  as  done 
previously. 

The  effect  of  the  cockpi t/radome  section  is  determined 
using  flat  plates  attached  to  the  fuselage  just  as  the  wings. 

Since  computation  time  in  analyzing  a real  cockpit  radome 
section  is  so  great,  it  is  apparent  that  a simple  model  is  needed. 
Originally,  a finite  bent  plate  approximation  of  the  cockpit  and 
nose  section  was  considered  because  of  an  experimental  investigation 
of  the  radome.  This  resulted  in  a truncated  fuselage  model  as 
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discussed  in  Reference  [31].  However,  the  radiation  pattern  cal- 
culated using  a flat  plate  model  cotnoared  very  favorably  with  the 
bent  plate  model.  Hence,  the  flat  plate  model  was  adopted  not 
only  for  its  simplicity  in  analysis  but,  also,  for  its  reduced 
computer  running  time.  Notice  that  the  plate  is  mounted  at 
arbitrary  angles  relative  to  the  cylinder  to  simulate  the  cockpit/ 
radome  section  as  seen  in  Fig.  34  (a).  This  is  an  improvement 
over  the  previous  roll  plane  model  where  the  wings  are  restricted 
to  be  horizontal. 

Similar  considerations  were  made  in  studying  the  vertical 
stabilizer  which  resulted  in  the  selection  of  a bent  plate 
model.  Thus,  the  scattering  effects  of  the  cockpit  and  vertical 
stabilizer  can  be  taken  into  account  very  simply.  The  theoretical 
model  as  shown  in  Fig.  34  (a)  is  an  illustration  of  the  analytical 
simulation  of  an  actual  aircraft.  With  all  the  necessary  mechanisms 
being  completed,  the  radiation  and  scattering  by  the  fuselage,  wings, 
cockpit,  and  stabilizers  can  now  be  calculated.  Finally,  the 
volumetric  pattern  of  a fuselage  mounted  antenna  is  obtained 
by  summing  the  direct  field  from  the  source  and  the  reflected 
and  diffracted  field  from  the  various  scatters. 

Before  it  can  be  used  to  determine  the  volumetric  pattern 
fcr  an  airborne  antenna,  tnis  newly  developed  solution  must  be 
tested  and  verified.  The  elevation  pattern  of  an  axial  slot 
mounted  on  a prolate  spheroid,  thus,  was  calculated  and  is 
presented  in  Fig.  37.  The  comparison  between  the  measured  and 
calculated  results  is  very  persuasive.  This  particular  problem 
is  an  effective  test  in  that  it  illustrates  the  accuracy  for 
tre  limiting  case  of  a small  nearly  spherical  object  as  opposed 
to  the  large  cylindrical  aircraft  model. 

To  illustrate  the  validity  of  the  new  solution,  the  principal 
plane  radiation  patterns  of  antennas  mounted  on  the  KC135  aircraft 
tr.at  was  previously  computed  are  now  computed  with  the  improved 
solution.  With  the  cockpi t/radome  and  vertical  stabilizer  being 
taken  into  account,  the  new  computer  model,  consisting  of  an  80" 
by  3.75"  right  serr.i -el  1 ipse  and  an  8"  by  3.75"  left  semi-ellipse, 
is  shown  in  Fig.  38.  The  cockpi t/radome  section  is  simulated 
by  a flat  plate  and  the  vertical  stabilizer  is  approximated  by  a 
bent  plate  model  for  better  simulation.  The  elevation  plane 
patterns  for  a short  monopole  mounted  forward  and  over  the  wings 

Shewn  in  Figs.  39  (a)  and  (b),  respecti vely.  The  patterns 
f-n  a circumferential  KA-band  waveguide  are  shown  in  Figs.  40  (a) 
and  (d).  Tne  results  for  an  axial  KA-band  waveguid  are  depicted 
in  Figs.  41  (a)  and  (b).  It  is  observed  that  the  comparison 
between  the  calculated  and  measured  results  is  very  favorable. 
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Figure  38b.  Computer  simulated  model  for  the  cross-section 
fat  antenna  location)  of  a KC-135  aircraft 
(front  view). 


86 


Figure  38c.  Computer  simulated  model  for  a KC-135  aircraft  (top  view). 


Figure  39a.  Elevation  plane  pattern  for  a A/4  monopole 
mounted  forward  of  the  wings  on  a KC-135 
aircraft  (with  radome  and  vertical  stabilizer 
included). 
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Figure  39b.  Elevation  plane  pattern  for  a A/4  monopole 

mounted  above  the  wings  on  a KC-135  aircraft 
(with  radome  and  vertical  stabilizer  included). 
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Figure  40a.  Elevation  plane  pattern  for  a circumferential 
KA-band  waveguide  mounted  forward  of  the  wings 
on  a KC-135  aircraft  (wi th  radome  and  vertical 
stabilizer  included). 
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Figure  40b.  Elevation  plane  pattern  for  a circumferential 
KA-bar.d  waveguide  mounted  above  the  wings  on  a 
KC-135  aircraft  (with  radome  and  vertical 
stabilizer  included). 
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Figure  41a.  Elevation  plane  pattern  for  an  axial  KA-band 
waveguide  mounted  forward  of  the  wings  on  a 
KC-135  aircraft  (wi th  radome  and  vertical 
stabilizer  included). 
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Figure  41b.  Elevation  plane  pattern  for  an  axial  KA-band 
waveguide  mounted  above  the  wings  on  a KC-135 
aircraft  (with  radome  and  vertical  stabilizer 
included). 
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The  discrepancy  displayed  in  aft  sectors  in  Figs.  14  and  15  no 
longer  exist  in  Figs.  39  and  40.  For  axial  slot  cases,  the 
back  lobe  is  picked  up  as  ‘compared  with  previous  results.  Notice 
that  the  nulls  of  the  monopole  for  the  theory  and  measurement  do 
not  align  together  as  one  can  see  in  Figs.  39  (a)  and  (b).  This 
can  be  attributed  to  some  extent  to  the  misalignment  of  the  monopole 
and  surface  normal  in  the  experimental  work. 

Finally,  the  roll  plane  radiation  patterns  for  the  new  model 
are  calculated.  In  each  case,  the  pattern  is  only  slightly  improved 
over  previously  calculated  results  as  shown  in  Figs.  27  and  29.  This 
is  attributed  to  the  fact  that  the  surface  rays  contribution  from  the 
elevation  plane  model  is  negligible  in  the  plane  of  interest.  In  any 
event,  these  results  demonstrate  the  validity  and  applicability  of  the 
new  solution  in  determining  the  radiation  patterns  for  general  fuse- 
lage mounted  antennas. 

To  further  demonstrate  the  versatility  of  this  solution,  the 
radiation  pattern  for  an  infinitesimal  monopole  mounted  on  the  bottom 
of  the  fuselage  of  a Boeing  737  with  a dielectric  radome  is  calculated. 
The  model  for  this  case  is  shown  in  Fig.  42  where  the  radome  effect  is 
simulated  by  a flat  plate  mounted  vertically  upward.  This  model  is 
adopted  because  of  the  similar  consideration  as  discussed  previously. 
Even  though  double  diffractions  are  not  included  in  the  calculation, 
the  resultant  pattern  is  still  in  very  good  agreement  with  the  measured 
result  as  shown  in  Fig.  43.  This,  again,  illustrates  the  capability 
of  the  solution  in  predicting  the  radiation  pattern  of  fuselage  mounted 
aircraft  antennas. 


S"*o 
■ £ o 

^ fO  (VJ 

o w - 

tl  II  II 

a 13  ]d 


T^l 

<*  Jr 

Nl  & 

I ^ 

1? 


\ I ' 

■aL' 


4->  40 

ro  4— 

<0 

“O  S- 
0)  cj 

4->  Sw 
C -r- 
=1  <0 
o 

E r-  . 

ro^ 

a>  r*.  cd 

r—  -»-> 
Q CD  rO 
aCn 
O — Q. 

C CD 

O O *-> 
E ca  (tj 

^3-  HJ  »♦- 

4-  <o 

o 

rO  >> 
O)  xa 
t-  cn 
o o o 

<4-  r—  CD 
<D  i— 
r—  t/>  c— 
CU  3 d) 
T3  <4-  “O 

o o 

E 0)  £ 
-C 

T3  D cn 

QJ  c 

4- >  M—  •#— 

rO  O O 
r—  JD 

3 E 

E O C 

•r—  -4— > O 
CO  -*->  •!— 

O 4-> 

5-  -Q  O 

CD  CD 

OJ  U) 
3 -C 

a.4->  <d 

E CO 

o c o 
o o c 


/f 


0 = 0° 


180* 


Figure  43.  Elevation  plane  pattern  of  a x/4  monopole  mounted 
at  station  222  on  the  bottom  of  a Boeing  737  air- 
craft. ($  = 0°  at  the  left;  <f  = 180°  at  the 
right). 


96 


SECTION  IV 


VOLUMETRIC  PATTERNS  OF  AIRBORNE  ANTENNAS 


The  theoretical  approach  to  analyze  the  complete  volumetric 
radiation  pattern  of  fuselage  mounted  antennas  was  developed  in 
the  last  section.  This  theoretical  solution  was  based  on  the  two 
•principal  plane  analyses  presented  earlier  and  it  has  been  veri- 
fied by  applying  it  to  various  cases  previously  calculated  using 
principal  plane  analyses.  The  results  show  great  improvement 
over  previous  calculations  based  on  comparisons  with 
experimental  measurements.  Hence,  the  new  volumetric  solution 
is  now  applied  to  analyze  the  volumetric  patterns  of  fuselage 
mounted  airborne  antennas. 

The  aircraft  model  of  most  interest  in  our  study  is  the 
Seeing  737  aircraft  in  that  extensive  experimental  work  is  avail- 
able. The  volumetric  patterns  of  this  aircraft  model  is  examined 
in  this  section. 

The  models  used  to  simulate  a Boeing  737  aircraft  are  shown 
in  Fig.  44,  where  the  rol 1 plane  cross-section  of  the  fuselage 
is  approximated  by  a 65.86"  by  43.3"  elliptic  cylinder  for  an 
antenna  mounted  at  Station  220  on  the  top  of  the  fuselage.  The 
elevation  profile  is  modeled  by  a composite  elliptic  cylinder 
with  a = 48.72",  b = 308.56",  and  b*  = 1307.04".  Notice  that  the 
cocxpit  nose  section  and  vertical  stabilizer  are  approximated  by 
finite  flat  and  bent  plates  mounted  obliquely  on  the  fuselage 
wnich  results  in  a simple  model  for  our  aircraft  as  shown  in  Fig. 


In  Fig.  44(c),  the  top  view  of  our  model  is  shown  to 
illustrate  the  finite  three-dimensional  effect  of  an  aircraft. 
The  dotted  line  indicates  the  width  of  the  cross-sectional 
(roil  plane)  cylinder  used  in  the  calculations.  The  models 
of  che  wings,  cockpit  nose  section,  and  stabilizers  are  also  pre 
seated.  The  geometry  is  taken  directly  from  the  three  principal 
views  of  the  aircraft  scale  model.  As  presented  earlier,  the 
coordinate  systems  used  for  both  the  elevation  and  roll  plane 
geometries  are  depicted  in  Fig.  44.  , 

The  radiation  patterns  for  a A/4  monopole  mounted  at 
Station  220  above  the  cockpit  on  a Boeing  737  aircraft  are, 
tier.,  calculated  using  the  model  just  described  with  the 
analysis  presented  in  the  last  section.  The  three  principal 
plane  results  are  shown  in  Figs.  45  to  47  and  found  to  be  in 
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Figure  44a.  Computer  simulated  model  for  the  fuselage  profile 
of  a Boeing  737  aircraft  (side  view).  The  antenna 
is  located  at  station  220  on  top  of  the  fuselage. 
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Figure  44b.  Computer  simulated  model  for  the  cross  section 
(at  antenna  location)  of  a Boeing  737  aircraft 
(front  view).  The  antenna  is  located  at  station 
220  on  top  of  the  fuselage. 


99 


r i- — 


180° 


Figure  45.  Elevation  plane  pattern  of  a a/4  monopole  mounted 
at  station  220  on  top  of  a Eoeing  737  aircraft. 

( y = 0°  at  the  left;  4 = 180°  at  the  right). 
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Figure  46.  Roll  plane  pattern  of  a a/4  monopole  mounted  at 
station  220  on  top  of  a Boeing  737  aircraft. 

($  = 90°  at  the  left;  $ = 270°  at  the  right). 
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Figure  47.  Azimuth  plane  pattern  of  a A/4  monopole  mounted 
at  station  on  top  of  a Boeing  737  aircraft. 

(0  = 92°), 
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very  good  agreement  with  measurements . The  experimental  work 
was  performed  with  a great  deal  of  patience  and  precision  by 
the  technical  staff  at  NASA  (Hampton,  Va.),  using  a 1/llth 
scale  model  of  Boeing  737  aircraft. 

The  coordinate  system  used  for  the  experimental  measurements 
is  shown  in  Fig.  48,  in  which  the  z-axis  is  vertically  upward. 

In  order  to  calculate  a radiation  pattern  in  terms  of  this 
experimental  coordinate  system,  a transformation  of  coordinates 
is  necessary  so  that  the  corresponding  radiation  direction  in 
the  analytical  aircraft  reference  coordinate  system  can  be 
determined.  This  is  accomplished  by  considering  the  z-axis  of 
the  experimental  coordinate  system  as  a radial  vector  in  the 
reference  coordinate  system.  By  inputting  its  spherical  coordin- 
ates  (&Zref  = 90° » ^Zref  = 0°)  as  the  rotation  angles  ec  and  $ 

required  in  the  computer  program,  the  corresponding  radiation 
direction,  in  terms  of  the  reference  coordinate  system,  is 
obtained. 

Z 


Figure  48.  Illustration  of  the  coordinate  system  used  for 
experimental  measurements. 
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The  complete  volumetric  pattern  for  a A/4  monopole  located 
at  Station  220  on  top  of  a Boeing  737  aircraft  is  calculated. 

The  off-principal  plane  elevation  patterns  are  shown  in  Figs. 

49  to  52.  The  azimuth  plane  patterns  are  shown  in  Figs.  53  to 
60.  In  each  case,  the  calculated  results  compare  very  favorably 
with  the  measurement.  It  is  noted  that  the  measurement  results 
have  some  asymmetry  in  the  patterns.  This  could  be  attributed  to 
misalignment  of  the  monopole  with  respect  to  the  surface  normal 
or  the  movement  of  the  model  due  to  shifting  weight  during  the 
measurement. 
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Figure  49.  Elevation  plane  pattern  of  a a/4  monopole  mounted 
at  station  220  on  top  of  a Boeing  737  aircraft. 

($>  = 10°  at  the  left;  $ = 190°  at  the  right). 
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Figure  50.  Elevation  plane  pattern  of  a a/4  monopole  mounted 
at  station  220  on  top  of  a Boeing  737  aircraft. 

($>  = 20°  at  the  left;  <P  = 200°  at  the  right). 
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Figure  51.  Elevation  plane  pattern  of  a >./4  monopole  mounted 
at  station  220  on  top  of  a Boeing  737  aircraft. 

(9  = 30°  at  the  left;  9 = 210°  at  the  right). 
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Figure  55.  Azimuth  plane  pattern  of  a A/4  monopole  mounted 
at  station  220  on  top  of  a Boeing  737  aircraft. 
(6  = 70°) 
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Figure  56.  Azimuth  plane  pattern  of  a x/4  monopole  mounted 
at  station  220  on  top  of  a Boeing  737  aircraft, 
(e  = 80 J ) 
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Figure  57.  Azimuth  plane  pattern  of  a a/4  monopole  mounted 
at  station  220  on  top  of  a Boeing  737  aircraft, 
(e  = 90°) 
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SECTION  V 


SURFACI  CURRENT  ANO  CHARE, I DINSITY  INOUCfn  ON  AIRCRAI  T MODUS 


Based  on  the  GTD  analyses  for  the  airborne  antenna  problems,  this 
chapter  is  devoted  to  calculate  the  surface  current  and  charge  density 
induced  on  various  aircraft  models  using  the  reciprocity  theorem.  Re- 
call that  in  Chapter  I the  relation  between  the  antenna  radiation  pat- 
terns and  the  induced  surface  current  and  charge  density  has  been  dis- 
cussed. For  convenience,  these  relations  are  rewritten  as  follows: 
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(74) 


Using  Eq.  (74)  and  the  GTD  analyses,  numerical  results  of  induced  surface 
current  and  charge  density  are  presented  for  various  aircraft  models. 

Consider  a Boeing  747  illuminated  by  a plane  wave  with  broadside 
incidence.  The  geometry  of  the  aircraft  and  its  computer  model  are  shown 
in  Fig.  61  with  observation  point  located  on  the  top  of  the  fuselage 
denoted  by  letter  P in  the  figure.  Note  that  the  fuselaqe  axis  is 
assumed  to  be  the  z-axis.  The  magnitude  and  phase  of  induced  current 
and  charge  densities  for  two  different  polarizations  are  presented 
in  Figs.  62-65.  The  THETA  MODE  and  PHI  MODE  appearing  in  these  figures 
indicate  that  the  electric  field  intensities  associated  with  the  in- 
cident plane  wave  are  o-polarized  and  ^-polarized,  respectively. 

The  computed  model  of  a KC-135  is  shown  in  Fig.  66.  A plane  wave 
is  incident  at  o=45°  and  ^=45°.  The  surface  current  and  charge  densities 
induced  on  a KC-135  aircraft  are  shown  in  Figs.  67-70. 
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Finally,  a Boeing  737  aircraft  is  investigated.  The  geometry  and 
its  computer  model  are  shown  in  Fig.  71.  The  observation  point  is 
located  at  station  22u  on  top  of  the  fuselage.  The  coordinate  system 
used  for  the  calculation  is  illustrated  in  Fig.  72.  The  complete  volu- 
metric pattern  for  the  charge  density  induced  at  station  220  on  top  of 
a Boeing  737  is  shown  in  Fig.  73. 


119 


Figure  61b.  Computer  model  of  a 747  aircraft 
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PLANE  NAVE  INCIDENT  AT  PHI-  0.00  DEGREES 

THI-  90.00  DEGREES 


PHI  MODE 

Figure  65.  The  surface  current  and  charge  densities 
induced  on  a 747  aircraft. 
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Figure  68.  The  surface  current  and  charge  densities 
induced  on  a KC-135  aircraft. 
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PLANE  HAVE  INCIDENT  AT  PHI-  US. 00  DEGREES 

THI-  1*5-00  DEGREES 


Fiqure  71.  Computer  simulated  model  for  a Boeing  737  aircraft 
(top  view).  The  observation  point  is  located  at 
station  220  on  the  top  of  the  fuselage. 
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SECTION  VI 

SUMMARY  AND  CONCLUSIONS 


High  frequency  solutions  for  the  three-dimensional  surface  current 
and  charge  density  induced  on  the  aircraft  fuselage  for  arbitrarily- 
incident  plane  wave  has  been  the  object  of  this  research.  It  is  well 
known  that  the  induced  surface  current  and  charge  density  due  to  plane 
wave  incidence  are  related  to  the  radiation  pattern  due  to  infini tesmal 
tnonopole  and  slot  antennas.  Since  extensive  experimental  scale-model 
measurements  for  flush-mounted  antennas  on  aircraft  are  available  for 
comparison,  the  analyses  for  the  induced  surface  current  and  charge 
density  is  carried  out  in  terms  of  antenna  problems.  This  report  pre- 
sents an  accurate  and  efficient  solution  for  the  problem  using  the  geo- 
metrical theory  of  diffraction  (GTD). 

The  basic  approach  applied  here  is  to  break  up  the  aircraft  into 
its  simplest  structural  forms.  Due  to  the  high  frequency  requirements 
of  airborne  antennas,  these  structures  may  be  analyzed  using  ray  optics 
techniques  with  numerical  values  obtained  using  GTD.  Once  the  scatter- 
ing from  these  various  structures  is  found,  it  is  then  adapted  to  the 
aircraft  model  simply  by  adjusting  the  field  incident  on  the  structural 
scatterer.  The  only  limitation  of  the  GTD  solution  is  that  the  source 
and  various  scattering  centers  be  separated  by  at  least  a wavelength. 

In  some  cases,  even  this  requirement  can  be  relaxed. 

A theoretical  solution  was  developed  in  Chapter  IV  to  analyze 
complicated  three-dimensional  volumetric  radiation  patterns  for  fuse- 
lage mounted  airborne  antennas.  This  solution  utilizes  roll  and  ele- 
vation plane  model  analyses  developed  earlier.  The  procedure  to  com- 
bine these  two  solutions  into  a complete  solution  that  can  handle  the 
volumetric  pattern  is  based  on  a previous  study  of  antennas  mounted  on 
prolate  spheriods.  The  belt  concept  of  blending  these  two  solutions 
together  is  the  key  to  the  success  of  the  complete  solution.  The  use 
of  flat  or  bent  plates  to  approximate  an  aircraft  cockpit  nose  section 
and  vertical  stabilizer  is  also  new  and  quite  useful,  not  only  because 
of  its  simplicity  in  analysis  but  also  due  to  its  practicality  in  com- 
putation. 

A numerical  procedure  that  can  be  used  to  model  a practical  air- 
craft fuselage  has  also  been  presented.  The  three-dimensional  fuselage 
in  the  numerical  model  is  simulated  by  two  composite  elliptic  cylinders; 
one  approximates  the  aircraft  profile  and  the  other  its  cross-section. 
These  ellipses  are  numerically  obtained  through  a best-fit  ellipse 
routine.  The  description  of  the  wings,  cockpit,  and  stabilizers  are 
measured  directly  from  the  three  principal  views  of  the  aircraft. 
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To  demonstrate  the  validity  and  applicability  of  the  theoretical 
solutions,  the  radiation  patterns  in  the  principal  planes  for  antennas 
mounted  on  various  locations  on  a KC-135  aircraft  and  a Boeing  737  air- 
craft have  been  calculated  and  compared  with  measurements  taken  at  NASA 
(Hampton,  Va.).  The  results  show  good  agreement  between  the  analytical 
and  experimental  work  and  obvious  improvement  over  previous  analses. 

To  illustrate  the  versatility  of  the  new  solution,  the  radiation  patterns 
for  antennas  mounted  on  a space  shuttle  with  the  radome  taken  into  account 
have  been  calculated.  Again,  the  result  is  very  encouraging.  The  off- 
principal  plane  radiation  patterns  of  a A/4  monopole  mounted  above  the 
cockpit  of  a Boeing  737  aircraft  have  been  calculated.  The  patterns 
obtained  were  shown  to  compare  well  with  measurements  taken  at  NASA 
(Hampton,  Va.),  which  verify  the  assumptions  made  in  the  theoretical 
solutions. 

Finally,  some  numerical  results  for  the  surface  current  and  charge 
density  induced  on  various  aircraft  models  are  calculated  and  presented 
in  Section  V.  The  solutions  developed  in  this  research  provide  a useful, 
accurate,  economical,  and  efficient  means  for  determining  the  surface 
current  and  charge  density  induced  on  aircraft  fuselage  due  to  arbitrary 
plane  wave  incidence.  For  example,  the  program  developed  in  this  study 
has  now  been  delivered  to  NASA  (Hampton,  Va.),  and  it  typically  runs  a 
pattern  in  30  seconds  or  less  on  a CDC  6600  digital  computer. 

It  should  be  noted  that  in  the  present  study,  the  effect  of  the 
engines  is  negligible  and  they  have  been  ignored  in  the  model.  However, 
the  engines  appear  to  have  some  effect  on  the  radiation  patterns  when 
an  antenna  is  mounted  on  the  bottom  of  the  fuselage.  When  necessary, 
the  engine  effect  can  be  included  in  the  solution  by  using  a finite 
elliptic  cylinder  as  presented  in  Reference  [55]. 


APPENDIX  A 

COORDINATE  SYSTEMS  TRANSFORMATION 


Consider  that  a vector  R (x,  y,  z)  in  an  xyz  coordinate 
system  as  shown  in  Fig.  A-la.  A second  coordinate  system, 
(x'y'z1)  is  also  depicted,  which  is  formed  by  first  rotating 
the  x-axis  of  the  original  coordinate  system  through  an 
angle  of  4>0  about  z-axis  and  then,  the  z-axis  through  an  angle 
of  e0  about  the  rotated  x'-axis^  Recall  that  in  a vector  space, 
the  coordinates  x,  of  a vector  X representing  a point  in 
n-dimensional  space  canalso  be  regarded  as  the  coefficients  of 
the  unit  vectors  ej  if  X is  represented  as  a sum  of  multiples 
of  the  unit  vectors 


X = 


- x,  e,  ♦ 


+ Xn  en 

n n 


(A-l) 


These  unit  vectors  form  a vector  basis  which  span  the  n-dimensional 
vector  space.  In  the  three-dimensional  geometry,  where  n = 3, 
let  {x,  y,  z}  be  a basis  of  the  vector  space  under  consideration. 
Thus,  the  vector  R can  be  written  in  terms  of  these  unit  vectors 
as 


R = xx  + yy  + zz 


(A-2) 


Let  {x‘,  y',  z ' } be  another  basis  of  the  same  space  such  that 
the  vector  R is  expressed  as 


R = x 1 x 1 + y'y'  + z ' z ‘ (A-3) 

Then,  by  a linear  transformation,  the  coordinates  of  the  vector 
R relative  to  the  {x,  y,  z}  b$siSAcan_be  determined  from  the 
coordinates  relative  to  the  {x‘,  y',  z'}  basis  through  a matrix 
relation  as 


X = PX’ 


(A-4 ) 
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Figure  A-la.  Illustration  of  coordinate  system  rotation. 

(Rotate  $ angle  about  z axis  and  then  rotate 
e0  angle  ibout  x'  axis). 


Figure  A-lb.  Illustration  of  coordinate  system  rotation. 

(Rotate  $ angle  about  z axis  and  then  rotate 
e angle  about  y*  axis). 
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Where  X and  X'  are  column  vectors  given  by 


and  PJs  a 3x3  transformation  matrix  from  the  {x,  y,  z)  basis  to 
lx',  y',  z'}  basis.  By  multiplying  Eq.^A-6)  by  P_l,  the  inverse 
of  P,  the  coordinates  of  R”  in  {x‘,  y',  z'}  basis  can  be  found  in 
terms  of  the  coordinates  in  {x,  y,  z}  basis,  i.e., 

X'  = P_1X  (A-6) 


Eqs.  (A-4)  and  (A-6)  describe  the  relation  on  which  the  coordinates 
of  a vector  in  one  coordinate  system  can  be  determined  from  its 
coordinates  relative  to  another  coordinate  system. 

Since^the  unit  vectors  x' , y' , and  z'  are  related  to  the  unit 
vectors,  x,  P,  and  2 in  the  following  manner,  i.e., 

x'  = cos  <t>Q  x + sin  $o  y 

* a * 

y'  = -cos  oQ  sin  <t>Qx  + cos  oq  cos  4>Qy  + sin  eQz  (A- 7 ) 

z'  = sin  sin*  x -sin  cos  $ y + cos  e z 
O O 0 0 0 

the  transition  matrix  ?x,  which  is  defined  as  the  transpose  of 
the  above  matrix  of  coefficients,  can  be  obtained  as 
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where  subscript  x indicates  the  e angle  rotation  is  about  the 
rotated  x'-axis.  By  inverting  the  matrix  Px,  the  inverse  matrix 
Px"'  of  Px  can  be  easily  determined  and  is  given  by 
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Similarly,  if  a new  coordinate  system,  x'y'z'  system,  is  set  up 
by  first  rotating  the  x-axis  of  an  angle  of  $0  about  z-axis  and 
then  the  z-axis  an  angle  of  e0  about  the  y'-axis  in  a counter- 
clockwise sense  as  seen  in  Fig.  A-lb,  the  transformation  matrices 
P and  P can  be  determined  as 


and 
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cos  0O  . 

where  the  subscript  y indicates  the  rotation  of  0O  angle  is  about 
the  rotated  y'-axis  in  distinction  with  that  of  the  x'-axis  in 
Fig.  A-la. 
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APPENDIX  B 


GEOMETRICAL  CONFIGURATION  FOR  OFF-CENTER  LINE 
MOUNTED  ANTENNAS 


Consider  a source  which  is  located  at  some  angle  off  the  center 
line  of  an  aircraft  fuselage  as  shown  in  Fig.  B-l . This  source  can  be 
an  infinitesimal  monopole  or  arbitrarily  oriented  slot  element.  As 
discussed  in  Section  III,  the  surface  geometry  nearest  to  the  antenna 
has  the  most  significant  effect  on  the  radiation  patterns.  For  this 
case,  the  longitudinal  surface  curvature  in  the  plane  which  is  parallel 
to  both  the  Zref  axis  and  the  surface  normal  at  the  source  location, 
and  the  transverse  surface  curvation  in  the  other  plane  which  is 
orthogonal  to  the  longitudinal  plane,  as  shown  in  Fig.  B-l,  play  a 
dominant  role  in  predicting  the  radiation  patterns  of  airborne  antennas. 
However,  the  transverse  surface  profile  which  cuts  obliquely  across 
the  fuselage  is  not  easily  obtained  from  a scale  model  drawing  of  an 
aircraft  (usually,  only  the  elevation  profile  and  cross  section  are 
given).  Thus,  to  obtain  the  necessary  elliptic  cylinders  for  the 
volumetric  pattern  analysis,  the  aircraft  cross  section  which  cuts 
orthogonally  through  the  source  location  is  used  to  approximate  the 
required  transverse  surface  profile  cylinder.  This  approximation  is 
reasonable  since  the  transverse  surface  profile  of  an  actual  aircraft 
does  not  change  drastically  for  a small  angle  deviation.  For  the  other 
elliptic  cylinder  which  approximates  the  longitudinal  surface  profile 
at  the  source  location,  the  following  approach  can  be  employed. 

First,  let  us  approximate  the  elevation  profile  of  the  aircraft 
of  interest  by  a composite  elliptic  cylinder  as  done  previously 
in  Section  III.  Its  necessary  parameters  are  given  by  a0,  semi- 
major axis,  b0  and  b0' , semi-minor  axis  for  left  and  right  half 
elliptic  cylinder  respectively,  as  shown  in  Fig.  B-2.  Recall  that 
the  source  location  in  our  volumetric  analysis  is  defined  by 
(<t>s,zs),  as  seen  in  Figs.  B-2  and  B-3,  relative  to  the  aircraft 
fuselage  cross  section  where  the  source  is  located.  Since  the 
position  of  the  source  is  given  in  terms  of  an  elliptic  cylinder 
approximating  the  fuselage  cross  section  at  the  source  location, 
its  Cartesian  coordinates  can  be  easily  determined  as 
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Figure  B-2.  Front  view  of  the  roll  plane  (transverse)  cylinder 
showing  the  antenna  location  and  Its  surface  normal. 
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ROLL  PLANE 
CYLINDER 


X$*  * P$  COS 


V = PS  sin  ♦s 


(B-l ) 


where 


af  bf 

2 2 2 n 

cos  4>s  + af  sin  <t>$ 


The  quantities  ar  and  br  are  the  semimajor  and  semiminor  axes 
of  the  cross  sectional  elliptic  cylinder.  In  terms  of  the 
elliptic  cylinder  coordinate  system  as  discussed  in  Chapter  V, 
x$'  and  y$'  are  given  by 

xs‘  = af  cos  v$ 


(B-2) 


where 


ys*  = bf  sin  v. 


»s  - ^ tan  *s 


From  Fig.  B-3,  the  distance  xeo  between  points  P and  Q is 
found  via  the  ellipse  equation  as 


b;  Jbo  ‘ zs  if  zs  is  ne9at1ve 


b1-  Jbo  " zs  if  zs  1S  P°s1t1ve 
o 


(B-3) 


where  a0»  b0,  and  b^,  have  been  defined  earlier.  Since  the  origin  of 
the  cross  sectional  cylinder  is  defined  at  from  point  P along 
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line  PQ,  the  source  position  In  terms  of  the  reference  coordiante 
system  can  be  easily  determined  from  Figs. 


x 


s 


- a^)  + a^  cos  v$ 


(B-4) 


4 


ys  ■ bf  sin  vs 


Once  the  source  location  is  determined  in  the  reference 
coordinate  system,  the  actual  composite  elliptic  cylinder  to 
simulate  the  longitudinal  surface  profile  at  the  antenna 
location  is  given  as 


xs  + ys 


ae  = < 


V 


2 2 
xs  + *s 


if  zc  is  positive 


(B-5) 


if  z$  is  negative 


be 


As  described  earlier,  the  elliptic  cylinder  which  simulates  the 
fuselage  cross  section  at  the  source  location  is  used  to  arroximate 
the  transverse  surface  profile.  Consequently,  the  necessary  axes 
of  the  transverse  elliptic  cylinder  are  given  by 


Note  that  the  transverse  surface  profile  discussed  previously 
lies  in  a plane  which  is  parallel  to  the  surface  normal  at  the 
source  location.  Hence,  the  transverse  or  roll  plane  elliptic 
cylinder  needs  a tilt  such  that  the  cylinder  surface  normal 
aligns  with  the  original  surface  normal  as  seen  in  Fig.  B-4. 
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The  required  tilt  or  rotation  (er  $r  tilt)  is  given  by 

r 


er  tilt  =' 


-1  ae  zs 

tan  — - — if  zf  is  positive 

2 12 2 s 

beH  + *s 


(B-7) 


tan 


-1  ae  2s 


. ,2 1 2 ” 2 

be  K + ys 


if  z$  is  negative 


V tilt 


- 0 


with  respect  to  the  reference  coordinate  system.  It  is  noted  that 
the  rotation  described  here  is  about  yre^  axis  and  not  xrefr  axis. 


Notice  also  that  the  coordinate  system  of  the  elevation  plane 
(longitudinal)  model  is  rotated  an  angle  (<j>n)  which  is  the  angle 
between  the  surface  normal  and  xref  axis  as  shown  in  Fig.  B~2 
relative  to  the  reference  coordinate  system.  In  other  words,  the 
elevation  model  coordinate  system  shown  in  Fig.  B-7  is  obtained 
through  a rotation  Hilt)  the  reference  coordinate  system 

as  discussed  in  Appendix  I.  The  rotation  angles  are  given  by 


ee  tilt 


90° 


*e  tilt  ” ^n 


(B-8) 


where 


I 


tan 


sin  v$ 
cos  v$ 


with  respect  to  the  reference  coordinate  system.  The  geometry  of 
the  coordinates  is  shown  in  Fig.  B-5. 


With  the  necessary  elliptic  cylinders  being  determined,  the 
procedure  to  combine  these  two  models  to  simulate  the  actual 
aircraft  is  the  same  as  discussed  in  Section  III.  Thus,  the 
volumetric  patterns  for  off-center  fpselage  mounted  airborne 
antennas  can  be  determined. 
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